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The DC electrical conductivity of an ultrarelativistic QED plasma is studied in real time by 
implementing the dynamical renormalization group. The conductivity is obtained from the real- 
time dependence of a dissipative kernel closely related to the retarded photon polarization. Pinch 
singularities in the imaginary part of the polarization are manifest as secular terms that grow in 
time in the perturbative expansion of this kernel. The leading secular terms are studied explicitly 
and it is shown that they are insensitive to the anomalous damping of hard fermions as a result of 
a cancellation between self-energy and vertex corrections. The resummation of the secular terms 
via the dynamical renormalization group leads directly to a renormalization group equation in real 
time, which is the Boltzmann equation for the (gauge invariant) fermion distribution function. A 
direct correspondence between the perturbative expansion and the linearized Boltzmann equation is 
established, allowing a direct identification of the self-energy and vertex contributions to the collision 
term. We obtain a Fokker-Planck equation in momentum space that describes the dynamics of the 
departure from equilibrium to leading logarithmic order in the coupling. This equation determines 
that the transport time scale is given by tti = -^^j^^^ij^- The solution of the Fokker-Planck 
equation approaches asymptotically the steady-state solution as ~ g-t/(4.038... ttr) "Yhe steady- 
state solution leads to the conductivity a = ^iin^i/^) to leading logarithmic order. We discuss the 
contributions beyond leading logarithms as well as beyond the Boltzmann equation. The dynamical 
renormalization group provides a link between linear response in quantum field theory and kinetic 
theory. 
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I. INTRODUCTION 



Transport phenomena play a fundamental role in the dynamics of the formation and eyolution of an ultrarelatiyistic 
quark-gluon plasma as well as in electromagnetic plasmas in the early uniyerse. The viscosity of a quark-gluon plasma 
enters in a hydrodynamic description 1, 2] and energy losses depend in general on transport coefficients of the 
plasma . In early universe cosmology the electric conductivity plays a fundamental role in the formation and decay 
(diffusion) of primordial magnetic fields I'll . 

A reliable calculation of transport coefficients in QED and QCD requires a deeper understanding of screening by 
the medium in order to treat the infrared properties of gauge theories systematically. A study of transport coefficients 
for hot and dense QED and QCD plasmas in kinetic theory including screening corrections has been presented in 
Refs. 0,0. In these references transport coefficients are computed by solving the Boltzmann kinetic equation for the 
single particle distribution functions with a collision term that includes screening corrections. In an ultrarelativistic 
plasma the propagators for soft quasiparticles have to include nonperturbative resummations in terms of hard thermal 
loops (HTL) [3, Q • This resummation includes consistently Debye screening for the longitudinal component of the 
gauge field (Coulomb interaction) and dynamical screening by Landau damping for the transverse component 
Screening via the HTL resummation of the propagators for gauge fields renders the transport cross sections infrared 
finite 

The program of calculating transport coefficients from a kinetic Boltzmann equation with HTL corrections to the 
collisional cross sections has been pursued further, and a number of transport coefficients have been calculated to 
leading logarithmic order in the coupling constant or in the limit of the large number of flavors @. An alternative 
approach to transport from a more microscopic point of view is based on Kubo's linear response formulation ^lOj . In 
this formulation, transport coefficients are related to two-point correlation functions of composite operators in the 
long- wavelength, low-frequency limit pTLll2||. This approach provides a direct link to transport phenomena from the 
underlying microscopic quantum field theory, and formally a correspondence between Kubo's linear response and the 
Boltzmann equation has been established in the literature for some scalar field theories 0, 0| . 

The calculation of transport coefficients from a microscopic quantum field theory at finite temperature or density 
involves non-perturbative resummations |12|| of a select type of Feynman diagrams. As discussed in Ref. [T^ for 
scalar theories this resummation is necessary because in the limit of vanishing external momentum and frequency 
there emerges pinch singularities corresponding to the propagation of on-shell intermediate states over arbitrarily 
long times and distances. This propagation will be damped by collisional processes which result in a perturbative 
width in the propagators. Each order in the perturbative expansion which would be naively suppressed by another 
power of the coupling, introduces a small denominator that is controlled by the width, thus cancelling the powers of 
the coupling in the numerators. Therefore, all terms in the perturbative expansion from a select group of ladder-type 
diagrams contribute at the same order |12| | . The sum of ladder diagrams leads to an integral equation which is similar 
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to a Boltzmann equation for the single-particle distribution function in scalar theories 0, IT^ , a result that was 
anticipated in Ref. (15i | within the context of transport in electron-phonon systems. 

Whereas the Boltzmann equation is a very efficient method to extract transport coefficients without complicated 
resummations, few attempts had been made to understand transport phenomena and to extract transport coefficients 
directly from the underlying microscopic quantum field theory in gauge theories. Recently, Kubo's formulation and 
the sum of ladder diagrams with a width for the fcrmion propagator was carried out for the case of the color 
conductivity '16"| confirming the results obtained via the kinetic approach '17*1. This program was extended to study 
the electric conductivity in an ultrarelativistic QED plasma (18, ,19J. In Ref. 18] it was shown that the sum of ladder 
diagrams is equivalent to the steady-state form (i.e., no time derivatives) of the linearized Boltzmann equation obtained 
in Refs. |^ Q in the leading logarithmic approximation. This result was confirmed in Ref. where it was also 
shown that new diagrams must be included to fulfill the Ward identities, but that these diagrams do not contribute to 
leading logarithmic order. While the important contributions of Refs. |16L lia . Il9l| established the equivalence between 
the ladder resummation and the steady-state form of the Boltzmann equation, a direct relationship between linear 
response, the quantum field theoretical approach to resummation of pinch singularities in real time and the time 
dependent Boltzmann equation in gauge field theories was not yet available. 

There is a fundamental interest in transport phenomena in ultrarelativistic plasmas which warrants an understand- 
ing of the main physical phenomena from different points of view. There are systematic resummation schemes in 
quantum field theory which could provide an alternative to the kinetic description or lead to a systematic calculation 
of higher order effects that are not captured by the kinetic approach. Furthermore, a microscopic approach should 
allow a clear understanding of when a kinetic description is valid, and also provide a framework to study strongly 
out of equilibrium phenomena outside of the realm of validity of kinetic theory. A resummation program that begins 
from the equations of motion for correlation function is the Schwinger-Dyson approach which leads to a hierarchy 
of equations for higher order correlation functions. Suitable truncations of this hierarchy justified by the particular 
physical case would lead to a systematic calculation of transport coefficients. Such program was initiated in Ref. pol | 
for the case of the DC electrical conductivity in a high temperature QED plasma. 

An alternative program to study relaxation and transport is based on a real-time implementation of the renor- 
malization group 22, 23j. Just as in the usual rcnormalization group, this approach is based on a wide separation 
of scales. In real-time these scales are the microscopic and the transport time scales, which are widely separated 
in weakly coupled theories. The dynamical rcnormalization group equations determine the evolution of correlation 
functions and expectation values on the long time scales. In particular, the Boltzmann kinetic equation can be in- 
terpreted as a renormalization g roup equation for the single-particle distribution function where the rcnormalization 
group parameter is real time |23.f23| . In Refs. [2^12^ it was pointed out that the pinch singularities that are ubiqui- 
tous in finite temperature field theory and that require the resummation of the perturbative expansion are manifest 
in real time as secular terms, namely terms that grow in time in the perturbative expansion of real-time kernels. The 
resummation of these pinch singularities in the Fourier transform of the correlation functions in the limit of small 
momentum and frequency leads to integral equations [T^ , while directly in real time this resummation is achieved 
by the dynamical renormalization group j22l |23| | . 

The main point of the dynamical renormalization group (DRG) approach is that time acts as an infrared cutoff, the 
correlation functions do not have singularities at any finite time since singularities only arise in the infinite time limit. 
In transport phenomena, as discussed above, pinch singularities arise from the propagation of on-shell intermediate 
states over arbitrarily long time (or distances). In the same manner as the usual renormalization group resums the 
infrared behavior of correlation functions in critical phenomena, the dynamical renormalization group resums the 
long-time behavior of correlators or expectation values f22|. In the case of the single-particle distribution functions, 
the DRG equation has been shown to be the Boltzmann equation . 

In this article we begin a program to study transport phenomena directly from the underlying quantum field theory 
in real time. The main goal of the program is to provide an understanding of transport and relaxation implementing 
concepts and ideas from the renormalization group description (either in real time or frequency and momentum 
space). Such program will lead to a description of transport phenomena much in the same manner as in critical 
phenomena and in deep-inelastic scattering, where physics at widely different scales is studied via the renormalization 
group. An example of this interpretation has been recently provided in the identification of the Altarelli-Parisi-Lipatov 
equations, which describe the evolution of parton distribution functions in invariant momentum transfer Q^, as DRG 
(or Boltzmann) equations j2Jj. Thus the dynamical renormalization group links transport phenomena with critical 
phenomena. Furthermore, a description of transport phenomena directly from the underlying quantum field theory 
in real time lead to further understanding of the validity of the kinetic description as well as to deal with strongly out 
of equilibrium phenomena where a kinetic description may not be suitable or reliable. 

Goals and main results of this article. In this article we focus on studying the DC electrical conductivity in 
an ultrarelativistic QED plasma. We study transport phenomena directly from quantum field theory in real time by 
implementing a resummation of the perturbative expansion via the dynamical renormalization group. This program 
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allows to establish a direct relationship between linear response, the resummation of pinch singularities and the time- 
dependent Boltzmann equation. The solution of the linearized dynamical renormalization group equation leads to 
the transport coefficients, in this case the DC electric conductivity. Along the way this program also establishes 
several important aspects: (i) a direct relationship between pinch singularities in the perturbative expansion in 
linear response and their resummation via the Boltzmann equation in real time, (ii) a direct identification of the 
single-particle distribution function with the nonequilibrium expectation value of a gauge invariant bilinear operator, 
(iii) the relevant approximations that determine the validity of the kinetic description, (iv) a direct identification of 
the self-energy and vertex corrections in quantum field theory and the different terms in the linearized Boltzmann 
equation, highlighting that the transport time scale is a consequence of the cancellation of the anomalous damping rate 
between the self-energy and vertex corrections, and (v) a pathway to include corrections to the leading logarithmic 
approximation and to the Boltzmann equation, with a firm ground on quantum field theory. 
The main results of this article are summarized as follows. 

(i) We begin by studying the real-time dynamics of an initially prepared magnetic field fluctuation as an initial value 
problem in linear response. We relate the DC conductivity to the time integral of a kernel closely related to the 
retarded photon polarization. A perturbative evaluation of this kernel in real time reveals secular terms, namely 
terms that grow in time and invalidate the perturbative expansion. These secular terms are a manifestation of 
pinch singularities in the Fourier transform of this kernel. We highlight how these secular terms are manifest 
in the perturbative solution of the Boltzmann equation in real time. This is an important and revealing aspect 
of the study of quantum field theory in real time that establishes a direct link between perturbation theory in 
the quantum field theoretical approach, secular terms (pinch singularities) and perturbation theory at the level 
of the Boltzmann equation. This study also illuminates the resummation of secular terms performed by the 
Boltzmann equation. 

(ii) In early studies of the DC conductivity in QED plasmas it was recognized that a subtle cancellation between self- 
energy and vertex corrections makes the conductivity insensitive to the anomalous fermion damping rate l25l . 
This cancellation for a QED plasma, originally noted in Ref. has since been found in many other contexts |2d 
[2^ l28ll29| and was explicitly shown to occur in the ladder resummation approach to extract the conductivity [iM 
Il9l | . Here we study the cancellation between the self-energy and vertex corrections with HTL photon propagators 
for the DC conductivity and transport phenomena directly in real time. An important consequence of the 
perturbative approach in real time is that the secular terms directly indicate at which time scale perturbation 
theory breaks down. This time scale is associated with the transport time scale beyond which a 
resummation program like the DRG must be used. Our analysis in real time clearly indicates that the cancellation 
of the contribution from ultrasoft photon exchange between the self-energy and vertex corrections makes the 
transport time scale insensitive to the anomalous damping rate of hard fermions [2^ . Within the framework 
of transport phenomena this cancellation is at the heart of the distinction between the quasiparticle relaxation 
time scale and the transport time scale, which makes the calculation of the DC electrical conductivity much 
more subtle and complicated than that for the color conductivity in QCD. For color transport there is no such 
cancellation and the leading contribution to the color conductivity arises from the anomalous damping of the 
hard quarks 16] . 

(iii) From the study in linear response, we identify the quantum field theoretical equivalent of the gauge invariant 
single-particle distribution function whose equation of motion is related to the dissipative kernel that leads to 
the conductivity. We obtain the equation of motion for this distribution function in perturbation theory and 
show that its solution features secular terms at large times. We then implement the dynamical renormalization 
group to resume the leading secular terms. The dynamical renormalization group equation for the single-particle 
distribution function is recognized as the time-dependent Boltzmann equation. Furthermore, we show that the 
anomalous damping of hard fermions is manifest as anomalous secular terms in the perturbative solution of 
the Boltzmann equation in the relaxation time approximation. The perturbative solution of the linearized time- 
dependent Boltzmann equation is shown to feature exactly the same linear secular terms as those found in the 
perturbative evaluation of the dissipative kernel. This study clarifies the cancellation of the anomalous damping 
between the different contributions in the linearized Boltzmann equation, and confirms the identification of the 
various terms in the linearized collision term with the self-energy and vertex corrections to the dissipative kernel 
in perturbation theory. This is yet another link between the perturbative framework in quantum field theory and 
that at the level of the Boltzmann equation, establishing directly the resummation implied by the Boltzmann 
equation. Our study clarifies directly in real time how the Boltzmann equation resums the pinch singularities 
found in perturbation theory. This detailed and clear link between quantum field theory and the Boltzmann 
equation cannot be extracted from the simplified equivalence between ladder diagrams and the steady-state 
Boltzmann equation and requires the time-dependent Boltzmann equation. 
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(iv) After recognizing that the leading logarithmic contribution to the collision term of the linearized Boltzmann 
equation is dominated by the kinematic region of momentum exchange eT ^ q T between particles with 
typical momenta p > T, we expand the collision kernel in powers of q/p to obtain a Fokker-Planck equation 
which describes the time and momentum evolution of the departure from equilibrium of the distribution function 
in the linearized approximation. The time dependent Fokker-Planck equation is solved by expanding in the 
eigenfuntions of a positive definite Hamiltonian. For late times, its solution approaches asymptotically the 
steady-state solution as ~ g-t/(4.038... ttr) ^ where ^tr = e''T^in(i/e) • solved analytically the steady-state Fokker- 
Planck equation for small and large momenta which describes the small- and large-momentum behavior of 
the departure from equilibrium in the steady state. We use these analytic asymptotic solutions to calculate 
numerically the DC conductivity. We find the leading logarithmic expression for the DC conductivity which 
agrees to within less than 0.1 % with the results of Ref. 0- Furthermore, the Fokker-Planck equation allows to 
establish contact with the variational formulation used in Ref. 9] . 

(v) We discuss the diagrams that are necessary to be included in the collision term to next to leading logarithmic 
order in the coupling, and the range of validity of the Boltzmann equation. It is pointed out that in order to 
go beyond the simple Boltzmann equation, terms that describe spin precession must be included in the set of 
kinetic equations. 

The article is organized as follows. In Sec. we introduce the real-time approach to extract the DC conductivity 
from the hydrodynamic relaxation of long-wavelength magnetic fields. The DC conductivity is determined by the 
time integral of a dissipative kernel directly related to the photon polarization. This kernel will provide the link 
between quantum field theory and kinetic theory. In Sec. IIIII we analyze the kernel that defines the conductivity 
to lowest order in the hard-thermal-loop approximation and describe the strategy to resum the perturbative series 
by extracting the leading secular terms in time. Then we study in detail the resummed one-loop self-energy and 
vertex in the hydrodynamic limit. In this section we make a connection with the results of Refs. pH 123. 123] for the 
real-time behavior of the fermion propagator and its anomalous damping in order to identify the secular terms in the 
dissipative kernel associated with the anomalous fermion damping. The Ward identity between the self-energy and 
vertex is shown to be fulfilled in this approximation. In Sec. II VI we compute the imaginary part of the resummed 
two-loop transverse photon polarization using the resummed self-energy and vertex obtained above. We extract the 
hydrodynamic poles (or pinch denominators) which in real time are manifest as secular terms that grow in time. Then 
we analyze the cancellation alluded to above and extract the leading logarithmic behavior in the gauge coupling of the 
leading secular terms. We discuss in detail which diagrams and which region of the exchanged momentum contribute 
to the leading secular term to leading logarithmic accuracy. 

In Sec.0we establish contact with the Boltzmann equation approach by identifying the single-particle distribution 
functions and obtain their equations of motion in perturbation theory. The perturbative solutions to these equations 
feature secular terms. The dynamical renormalization group is introduced to resum the perturbative equation of 
motion, and the DRG equation is identified with the Boltzmann equation. We establish direct contact between the 
linearized Boltzmann equation and the results obtained from the perturbative quantum field theory of the previous 
section. We obtain a Fokker-Planck equation, we solve it in an eigenfunction expansion and find its asymptotic late 
time behavior yielding a steady-state solution. ^From it we extract the leading logarithmic contribution to the DC 
conductivity. In this section we also discuss the regime of validity of the Boltzmann approach and the contributions 
that must be included to go beyond leading logarithmic order as well as beyond the Boltzmann equation. Finally, we 
present our conclusions and a discussion of future avenues in Sec. IVll Two appendices are devoted to summarizing 
the imaginary-time and real-time propagators used in the main text. 

II. LINEAR RESPONSE 

In this section we relate the conductivity to the relaxation of the gauge mean field as well as to the current induced 
by an external electric field using linear response. The reason for delving on linear response is to recognize the 
main real-time quantity that leads to the conductivity and is the link between Kubo's linear response, the dynamical 
renormalization group and the Boltzmann approach. 

The Lagrangian density of QED, the theory under consideration, is given by 

C^^{i0-e4)i^-^F,,F'^-' , (2.1) 

where the zero-temperature mass of the fermion m has been neglected in the high temperature limit T ^ m. We 
begin by casting our study directly in a manifestly gauge invariant form (see also Ref. |22|). In the Abelian case it is 
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straightforward to reduce the Hilbert space to the gauge invariant states and to define gauge invariant fields. This is 
best achieved within the canonical Hamiltonian formulation in terms of primary and secondary class constraints. In 
the Abelian case there are two first class constraints: 

7ro = 0, V • TT = -eV''''V' , (2.2) 

where ttq and tt = — E are the canonical momenta conjugate to A° and A, respectively. Physical states are those 
which are simultaneously annihilated by the first class constraints and physical operators commute with the first 
class constraints. Writing the gauge field in terms of transverse and longitudinal components as A = A^ + At with 
V X Al = V • At = and defining 

*(a;) = ?A(x) e'^/ '^'2'"^-G(x-y).Az,(y) ^ (2.3) 

where G(x — y) the Coulomb Green's function satisfying V^G(x — y) = 6^{x — y), after some algebra using the 
canonical commutation relations one finds that AT(a;) and ^{x) are gauge invariant field operators. 

The Hamiltonian can now be written solely in terms of these gauge invariant operators and when acting on gauge 
invariant states the resulting Hamiltonian is equivalent to that obtained in Coulomb gauge. However we emphasize 
that we have not fixed any gauge, this treatment, originally introduced by Dirac is manifestly gauge invariant. The 
instantaneous Coulomb interaction can be traded for a gauge invariant Lagrange multiplier field which we call 
leading to the following Lagrangian density p3| 

£ §(i ^ - e-f°A° + 67 • At)* + ^ [(9^ At)^ + {VAy] . (2.4) 

We emphasize that should not be confused with the temporal gauge field component. 

The main reason to introduce the gauge invariant formulation is that we will establish contact with the Boltzmann 
equation for the single particle distribution function which must be defined in a gauge invariant manner. 



A. Relaxation of gauge mean field 



The strategy is first to prepare the system at equilibrium in the remote past and then to introduce an adiabatic 
external source Jcxt coupled to the transverse gauge field At. The external source will induce an expectation value 
for the gauge field, representing a small departure from equilibrium. The external field is switched-off at f = and 
the induced expectation value relaxes then towards equilibrium. The real-time dynamics of relaxation is studied as 
an initial value problem and the conductivity is extracted from the relaxation rate of long-wavelength perturbations. 
This approach has already been used in Refs. ,22, 23], where more details can be found. 

Introducing an external source Joxt with the following time dependence 

J:,t(x,i) = J:.t(x)e^*e(-0, (2.5) 
we find the equation of motion for the expectation value aT(x, t) = (At(x, i)) j^^^ to be given by 

a^(x,t)- V2a^(x,i)-H J dVn^^,t(x-x',i-0 4(x',t') = '/cxtWe'*e(-i) , (2-6) 
where Hj'jj(x — x', i — <') is the retarded polarization |2^ 

n;it(x-x',t-<') = -z([j^(x,t),j^(x',0])e(i-t') 

/ + CXD 7 
— ImW^ {uj, k) e-«[-'(*-t')-k-(x-x')] _ i') ^ (2.7) 
-oo 



= I 



where J — e*7* is the electromagnetic current and fc = |k|. Using the Fourier representation of 8(t — i'), we find 
that the spatial and temporal Fourier transform of the retarded polarization can be written in the form 

H-^(c.,fc)^P^^""'^(^"'^). (2.8) 

Taking the spatial Fourier transform of the equation of motion l|2.6|) and denoting the spatial Fourier transforms of 
the expectation value of the gauge field and the external source term as aT(k, <) and Joxt(k, i), respectively, we find 
that for t < the solution of the equation of motion (|2.t)|l with the source eq. H2.5|l is given by 

aT(k,i < 0) = aT(k,0) e"* , (2.9) 
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where aT(k, 0) and Jcxt(k) are related by eq. (|2.6|l for t < 0. In the limit e ^ this solution entails that aT(k, t ^ 
0) = 0. Introducing the function, 

g'^{k,t) = - lmW^{ko,k) e-*^-°* , (2.10) 

carrying out an integration by parts in the last term of the equation of motion (|2.6() and changing variables t — i' — > i' 
we find that for t > the equation of motion H2.6|l becomes 



i^(k,t) + [k^s'^ + g'^{k,o)] a^{k,t) ~ f dt' g'^{k,t') 4(k,i-i') 



0, (2.11) 



where we have used aT(k, i < 0) = in the limit e ^ 0. Eq. (|2.11|l manifestly describes the dynamics of the induced 
expectation value as an initial value problem in real time, which can be solved via Laplace transform once the kernel 
g^^{k,t) is determined. An important aspect of this equation is that it illuminates the connection with relaxation and 
dissipation. 

// the kernel t/*-'(fc,t) is localized in time in the region Q < t < imcm(fc), where imom(fc) determines the memory 
of the kernel, then for t ^ imcm(fc) we can expand in derivatives aT(k, t — t') — sltO^, t) — sltO^, t) t' + ■ ■ ■ and the 
equation of motion H2.6|l becomes an infinite series of higher time derivatives, 

a^(k, t) + [k^P^ + g'^k, 0)] a^{k, t) + a'^{k) ^(k, t) + (k) di^{k, t) + ■ ■ ■ ^ . (2.12) 

In the above expression we have introduced, 

a'^{k) ^- dt g'^{k,t) , cri^(fc) = / dt' t' g'^{k,t') , (2.13) 
Jo Jq 

where the upper limit in the time integrals is taken to infinity since by assumption t 3> imcm(fc)- 

We will now focus on the relaxation of gauge mean field in the long- wavelength limit. For fc — s- we expect that 
C/'-' (0,t) = g{t) (5*-' , with g{t) = ^J2i=iS^^iO,t). Hence, we define the DC electrical conductivity as 

-1 3 poo 

'^ = ^1™E'^^^(^) = - / dt'g{t'). (2.14) 

The function C/*^ (fc, 0) has a simple interpretation g^^ (fc, 0) — Ren*-' (0, fc),^ which can be understood from the spectral 
representation of the Fourier transform of the polarization eq. (|2.8|l . In the Matsubara representation, the inverse 
transverse photon propagator in the static limit {un = 0) is given by fc^ + IlT{k,0). Thus the fact that there is 
no magnetic mass in thermal QED |S| entails that limfc^g R-enT(0, fc) = 0. As a result we can neglect t/*-' (fc,0) in 
eq. (|2.12|l since only its transverse component enters the equation of motion for a^. Therefore, eq. (|2.12l) becomes 

aT(k, t) + a slt {k,t) + k^ ariKt) = , (2.15) 

with the solution for the long-time relaxation of small departures from equilibrium for transverse gauge fields of long 
wavelengths fc <C cr, 

aT(k,t) - 6-'='*/" . (2.16) 

This purely diffusive relaxation is the hallmark of slow decay of magnetic fields in the hydrodynamic limit, an important 
aspect of the generation and evolution of cosmic magnetic fields 4J. 

The validity of the derivative expansion leading to the local equation of motion (|2.15l) can now be assessed. Since 
slt = — (fc^/c)^ ar, the term ar(k, t) in eq. H2.15|l is subleading for k <^ a. Furthermore, the coefficient of the second 
derivative term in the derivative expansion eq. H2.12|l is given by 

dttg'^{k,t)^a'^{k)t„,,„,{k) . (2.17) 







^ Since lmTlx{ko, k) is an odd function of fco, it vanishes as —Of^ko for fco — * and there is no need to append a principal part prescription. 
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Hence, the second derivative term in the derivative expansion will be much smaller than the first derivative term 
displayed in eq. 12.12|l when 

« 1 . (2.18) 

a 

Anticipating that after a resummation program and up to logarithms of the electromagnetic coupling a cx T / a and 
that imcm is the transport relaxation time (to be confirmed later) with tmcm oc l/a^T, the validity of the derivative 
expansion is warranted for long wavelength k <C eT. In this regime we can also drop the usual kinetic (second 
derivative) term since k <^ a and the long-time limit is completely determined by the hydrodynamic form eq. H2.15|l . 
The conductivity a is also identified with the inverse of the magnetic diffusion coefficient, which in turn determines 
the Reynolds number in the equations of magnetohydrodynamics and the decay of magnetic fields in cosmology. 
Introducing a convergence factor e ^ 0+ in the definition of crfc in eq. (|2.13|) , we find 

a^k) --Td^ S^^i^^t) . . r I?^^ . (2.19) 



(2.20) 







TT fco(fco-ie) 

Using the dispersive representation eq. H2.8|l for the retarded polarization, we find that 



j=0 



a^^k) = - [n'^'(w,fc) - ffJ'(k,0)] 
Since Imn(ci;, k) and Ren(aj, k) are odd and even in oj, respectively, we find, 

' ^= 3|ToE^''(^)- (2.21) 



^,(^^^„imn:W) 



Eq. (|2.21|l is the well known result of Kubo's linear response, and the main point of revisiting it here is the relationship 
between the usual result and the time integral of the kernel Q which will be shown to have a simple correspondence 
with the Boltzmann approach. 

B. Induced current 

Consider introducing an external electric field £ = —Jix with switched-on adiabatically 

A'{^, t) = A'{x) e'* e{-t) for t < , (2.22) 

in the Lagrangian density this corresponds to the shift — > Ay + >A.t- In linear response, the induced current is 
given by 

(JXx,i)) =i J d^x' ([J*(x,t), J^'(x',t')]) Mri^',t') Oit-t') . (2.23) 

In terms of the polarization eq. H2.7|) . integrating by parts in time and taking the spatial Fourier transform on both 
sides we find 

{r{k,t))^~ f dt' g'^{k,t-t')£^{k,t') . (2.24) 

Since the background gauge field At is switched-on adiabatically, the external electric field £ vanishes for t < 0. 
Furthermore assuming that the external electric field is constant in time for t > 0, we find 

(J'(fc,i)) = - / dt' g''{k,t') £^{k) , (2.25) 
Jo 

where we have relabelled t' t — t' . At this stage it is convenient to define 

jHt) = {j\k = Q,t)), g'^t) = g'^k = o,t) , r = r(fc = o), (2.26) 
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hence for a spatially constant external electric field the induced drift current at asymptotically long times is given by 

r{t^oo)^<T£\ (2.27) 
where we have used eq. H2.14|l . The important aspect of eq. (|2.25(l is that 

±r(^t) = -g^^t)£\ (2.28) 

a relation that will allow us to establish contact with the Boltzmann equation (see Sec. 0. In the limit fc ^ 0, it 
must be that 

g^^t) = g{t) , git) = 1 g^^t) , (2.29) 

in terms of which one obtains 

g{t) dt . (2.30) 

The linear response analysis both for the relaxation of a gauge mean field as well as for the induced current clearly 
suggests that the important quantity to study in real time is the long time behavior of the kernel g^-'{k,t) which in 
fact will be the link with the kinetic approach to be studied later. 



C. Strategy to extract DC conductivity 

Eq. I|2.21|l is the usual definition of the conductivity via Kubo's linear response The computation of the 

conductivity extracted from the long-wavelength and low-frequency limit of the imaginary part of the polarization 
features infrared divergences in perturbation theory. These divergences are manifest as pinch singularities and are a 
consequence of the propagation of intermediate states nearly on-shell for arbitrarily long time and distance T^l . In- 
cluding a width for the particles via a partial resummation of the perturbative series regulates these pinch singularities 
but each power of the coupling constant is compensated by powers of the (perturbative width) in the denominator 
implying the non-pcrturbative nature of the transport coefficient and requiring a resummation scheme |l2l Il6l Il8| . 

Instead of following this route, we propose here a different approach. The discussion above highlights that in the 
real-time framework to extract the relaxation of long-wavelength fields the important quantity is the kernel g{k,t) 
defined by eq. H2.10|l . In particular the DC conductivity eq. (|2.14|) will be finite if the asymptotic long time behavior 
of g{k,t) is such that the total time integral is finite for fc ^ 0. This requires that the kernel g{k,t) be localized 
in time, namely it has short memory in the long-wavelength limit. The important aspect of focusing on the kernel 
g{k,t) is that it is finite for any finite k and t. Infrared divergences can only appear in the infinite time limit which 
will result in an infinite integral in time. The finite time argument plays the role of a regulator, in fact for finite time, 
the intermediate states can only propagate during this finite duration of time and the pinch singularities are therefore 
regulated |2^|23. These singularities will emerge in the long-time limit in the form of secular terms, i.e., terms in the 
perturbative expansion that grow in time |2^ [23 ■ This situation is similar to that in critical phenomena and in field 
theory, the perturbative series for the n-point functions is finite for a finite momentum cutoff, only when the cutoff is 
taken to be very large does the perturbative expansion diverges. 

The dynamical renormalization group introduced in Refs. |2^|2^ provides a resummation of the secular terms and 
leads to an improved asymptotic long-time behavior of real time quantities (see Refs. IS^,!!^ for details and examples). 
Thus the strategy that we propose is the following: (i) We first obtain the perturbative expansion of Imn*^ (a;, k)/Lo 
for Lj, fc — + 0. This perturbative expansion will feature singular denominators of the form (k • p — cj)^", where p is 
some loop momentum to be integrated out and n > 1. The Fourier transform in time to obtain g^-'{k,t) of these 
denominators will feature terms of the form t"^^ which are secular for n = 2, 3, . . . , signalling the breakdown of the 
perturbative expansion at long times. Thus the kernel g{k, t) is finite for any finite time and the pinch singularities are 
manifest in the infinite time limit, (ii) After extracting the leading secular terms (terms that grow the fastest in time 
at a given order in perturbation theory), we will use the dynamical renormalization group resummation introduced 
m Refs. '22', "2^ to resum these secular terms and improve the asymptotic long-time behavior. After this real time 
renormalization group resummation the kernel g{k,t) will have an improved and bound long-time behavior and the 
conductivity can now be extracted via the time integral of the kernel in the zero-momentum limit. This is akin to the 
resummation and improvement of the perturbative series provided by the usual renormalization group in field theory 
and critical phenomena. 
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As we will see in detail below, the resummation of the secular terms via the dynamical renormalization group 
provides the bridge between linear response and the Boltzmann equation in real time. We now carry out this program, 
first by extracting the leading secular terms in perturbation theory and then invoking the dynamical renormalization 
group resummation, which will lead to the Boltzmann equation. 

III. PERTURBATION THEORY 

This section is devoted to analyzing the HTL contribution to the photon polarization and the one-loop hard fermion 
self-energy, hard fermion-soft photon vertex with HTL-rcsummed soft internal propagator as well as the corresponding 
Ward identity in the hydrodynamic limit. 

A. Polarization in the hard thermal loop approximation 

While there are several alternative methods to obtain the imaginary part of the polarization, we will carry out our 
perturbative study in the imaginary-time (Matsubara) formulation of finite temperature field theory The most 
transparent manner to compute higher order corrections is to introduce dispersive representations for propagators 
and self-energies. To one-loop order, the leading order in the high temperature limit (hard thermal loop) contribution 
is obtained by using the free field fermion propagators in the one-loop polarization. The polarization is found to be 
given by 



W={ruu,k) = e'Tj2T^T I ^ ^ S{P) 7^ S{P + K) , (3.1) 



where K = (wfc,k) and P = {iujm,p) with and ujrn being the bosonic and fermionic Matsubara frequencies, 
respectively (see Appendix IXI for notations) . Performing the Matsubara frequency sums (see Appendix IXj) and then 
taking the analytic continuation ivk lu + iO, we find in the HTL limit k,LU <ti T 

Imn™^*^ (c^, fc) = 27Te^iu [ f jP [<5(k . p - ,,) J(k • p + c.)] , (3.2) 

J {liry dp 

where in the product of fermionic spectral functions only the terms of the form p+ and p_ p_ [see eq. (jA2ll ] 
contribute in the limit fc,cij — > 0. Hence, we find in the HTL approximation 



HTL.,(,) ^ 4 ^2 / f ^ ^ ^S.^ . (3.3) 



(27r)3 dp P P - '^^1 9 
^From the linear response relation eq. H2.28|l . we find in the hard thermal loop approximation 



r d^n p'^T'^ 
J\t) = -46^ j J2^P'^ • Vpn^(p) = ^£'- (3.4) 



It is clear that to this order the induced drift current will grow linearly in time. This is a result of the fact that this 
lowest order approximation does not include collisions. As we will see later, this low order result will emerge also in 
the Boltzmann equation (see Sec. Iv)! . 

B. Resummed one-loop self-energy of hard fermions 

For massless fermions, it proves convenient to decompose the self-energy into positive and negative helicity compo- 
nents as 

S(po,p) =7-(p) S+(po,p)-l-7+(p) S_(po,p) , I]±(po,p) = ^Tr[7±(p)E(po,p)], (3.5) 

where 7±(p) = (7° =F 7 • p)/2 are the positive and negative helicity projectors, respectively. The resummed one-loop 
self-energy S(P) of a hard fermion with momentum p ~ T, as depicted in Fig.^ it receives contributions arising from 
soft internal photon and soft internal fermion lines. Let us write 

E(P) = I]^P(P) + E^f(P) , (3.6) 
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(b) 



FIG. 1: Resummed one-loop hard fermion self-energy with (a) soft internal photon and (b) soft internal fermion. The line with 
a dot denotes the HTL-resummed propagator. 

where the superscripts 'sp' and 'sf denote the soft-photon and soft-ferniion contributions, respectively. Using the 
imaginary time propagators given in Appendix ^ we find 



where P — {iuJm,p) with p ^ T , Q — (w„,q), Q' = (iw„,q) with q -^T, and V^iti) = 5^-' — q^q-' is the transverse 
projector. For later convenience, we have assigned the loop momentum Q to the soft internal propagators. Fur- 
thermore, we will neglect the instantaneous Coulomb interaction here and henceforth since it does not contribute to 
the imaginary part. After performing the Matsubara frequency sums, we can rewrite the self-energy as a spectral 
representation 



This spectral representation will be useful to carry out the Matsubara sums when the self-energy is inserted in the 
polarization. Decomposing ImS(po,p) onto 7±(p) analogous to eq. ^A.b\ and using the properties (see Appendix 
A for the spectral densities) *PL,T{-qo,q) = -*PL,T[qo,q), *P±{-qo,q) = *Pt('?o,9), riB{-qo) = -[l + "s(9o)] and 
npi—qo) = 1 ^ np{qo), we obtain 





(3.9) 



ImE±(po,p) = Imi;!^(po,p) + Imi;!^(po,p) 



(3.10) 
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where 

ImS!^(po,p) = Tre^y ^2;;^ J dqo [uBiqo) + npUp + q\)] 

x{[Ki{p,q) *PL{qo,q) + K^{p,(i) *PT{qo,q)] <5(po + go - |p + q|) 

+ [K^ip,q) *PLiqo,q) + K^{p,ci) *PTiqo,q)] Sipo - qo + \p + n])} , (3.11) 
T ^sf( ^ 2 f d\ l + nB(|p-q|)-nf(go) 

x[^^(P,q) P+(go,g) +i^^(p,q) *P-iqo,q)] S{pQ - qo - |p + q|) 

+ [X^(p,q) >+(go,g)+X±(p,q) >-(go,-z)] <5(po + 9o + |p + q|)} , (3.12) 



with 

i^f(P,q) = i[l±p-?T^], if2±(p,q) = lT(p-q)(pT^-q), ifg^lp, q) = 1 T (p • P^)(q • P^) (3.13) 

^From the above expressions, one finds that Iml]±(— po,p) = Im^T (Po ? p) • If IniS(po,p) is finite on the particle 
and antiparticle mass shells po = ±p, respectively, then T{p) = ImS±(±p,p) is the damping rate for the hard 
(anti)fermion Q. 

C. Anomalous damping of fermions 

The leading contribution to the imaginary part of the hard fermion self-energy arises from the ultrasoft region of 
the transverse photon exchange for which q <^ eT [21L 123 . |2^ . This is so because in a high temperature QED plasma 
the fermionic excitation receives an effective thermal mass to/ = eT/VS and the instantaneous Coulomb interaction 
is Debye screened by the electric (or Debye) mass ujd — eT/VS, whereas the magnetic interaction is only dynamically 
screened by Landau damping. Hence, we will focus only on the soft-photon contribution ImS!^ in this subsection. 

In the ultrasoft region of the loop momentum q, q^ <C eT in eq. H3.11|l . we can approximate the distribution 
"-5(90) + it-fHp + q|) — T/qo and the HTL spectral function for the transverse photon |2lLl23.l23| 

*_PT{qo^ ^ 6iq^ 

qo q^ 

which as pointed out in Ref. "^l*! can be interpreted as the exchange of a magnetostatic transverse photon. Furthermore, 
in this ultrasoft limit and for hard fermion momentum p > T we can replace p -|- q ~ p. After some algebra, we find 
near the particle and antiparticle mass shells po ~ ±p |2l|, l23l | 

ImS±P(po w ±p,p) ^aT dq dx {1 - x^) <5(po T P - qx) , (3.15) 

Jo J-i 

where x — p The integration is straightforward and yields the following leading order result 



ImE'±P(39o « ±P,p) = aTln 



PoTP 



0(boTp|°) . (3.16) 



While the fermion damping rate is ill-defined, several resummation schemes, either based on the thermal eikonal 
(Bloch-Nordsieck) approximation that sums rainbow diagrams |2l| or via the dynamical renormalization group j23 . l23l |. 
lead to the conclusion that the fermion propagator at large times decays in real time as 

(*(x, t)*(0, 0)) cx e"" ^ * , (3.17) 

which corresponds to an inverse relaxation time scale for the hard fermion 

r = arin(l/e). (3.18) 

This is known as the anomalous fermion damping 25] and is determined by the kinematic region of ultrasoft transverse 
photon exchange with momentum q <ti eT. The region q > eT of transverse photon exchange leads to a subleading 
contribution to the damping of hard fermions of the order a^T [I8LIT9I I2H. 
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(b) 

FIG. 2: Resummed one-loop hard fermion-soft photon vertex with (a) soft internal photon and (b) soft internal fermion. The 
line (vertex) with a dot denotes the HTL-resummed propagator (vertex). 

The reasons to repeat this analysis here are twofold: (i) to clarify in real time that the anomalous damping rate 
eq. ()3.18(l does not contribute to the conductivity, as a consequence of a cancellation between self-energy and vertex 
corrections as anticipated in Ref. [25j | , and (ii) we will study the solution of the Boltzmann equation in the relaxation 
time approximation which will also feature the anomalous time dependence. 

As it will be discussed in detail below in Sec. lIVDl there is a precise cancellation of the contribution from ultrasoft 
photon momentum exchange between the self-energy and the vertex corrections to the polarization in perturbation 
theory. This cancellation will also be manifest in the linearized Boltzmann equation, and the analysis presented above 
will lead us to a direct identification of self-energy and vertex contributions in the linearized Boltzmann equation. 

D. Resummed one-loop fermion-photon vertex 

We now evaluate the resummed one-loop hard fermion-soft photon vertex, which consistently with the Ward 
identity also receives contributions from soft internal photon and soft internal fermion (see Fig. Writing, 



r''(p + K,p) = pp^''(p + K,p) + r^'f'iP + K,p) , 



(3.19) 
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we find in the imaginary time formalism, 

er^P'^'{P + K,P) = e^'Tj^ j s{P + Q)-i^ S{P + K + Q)f*DL{Q) + 

+ y S{P + Q) r S{P + K + Q) 7^ *Dt{Q) VSici)] , (3.20) 

+ Y*S{Q')*T'^{Q' + K,Q')*S{Q' + K)j^ Dt{P-Q)V^{^^)] , (3.21) 

where *r'^ is the HTL-resummed fermion-photon vertex, K = (ivk,^) with k <^T, and the rest of the four-momenta 
are the same as those defined for the self-energy. 

Indeed, using the Ward identities satisfied by the respective free and HTL-resummed quantities 

K^,r - S-\Q) - S-\K + Q) , *T^'[Q' + K, Q') = *S-\Q') - *S-\Q' + K) , (3.22) 

and the expressions for the vertices given by eqs. (|3.2U|) and H3.21|l . one can easily show that the resummed one-loop 
self-energy and vertex satisfy the Ward identities 

Kf, pp^*'(p + K,p) = j:'p{p + k)- j:'p{p) , r^^^iP + k,p) = j:'^{p + k)- t!'\p) , (3.23) 

hence 

if^r^(P + X,P) = S(P + if)-E(P). (3.24) 

We note that these Ward identities for the resummed one-loop self-energy and vertex are exact in the kinematic region 
under consideration p ^ T and k <ti T. 

Let us first concentrate on the soft-photon contribution F'^P'''. The sum over the Matsubara frequencies can be 
done straightforwardly. Using the identity 



ab \b — a J \a b 

where the principal value (PV) prescription is necessary to define the limit 6 — a — > 0, we find that, 

f dPa 

r^P-^(P + K, P) ^ J J dqo j dpo j dso [7' Pf{Po, P + q) 7^ Pf{so, P + k + q) 7 

X *PT{qo,q) "Prlq) +7° PF(po,P + q) 7^ pf(so, P + k + q)7° *pL((?o,g)] 
"s('7o) + "-f(po) nBiqa) +nF{so) 



So -Po - i^k 



PO - qo - So ~ qo ~ «W„i - Zl/fc 



(3.26) 



where the denominator l/(so — po — ivk) in the above expression should be understood with the principal value 
prescription as in eq. (|3.25|) . 

The kinematic region of interest for the polarization and the conductivity corresponds to hard external fermion 
p ^ T and soft external photon k <^T . Furthermore after the analytic continuation of the external photon frequency 
iuk — > w -|- «0 we need u> -^T, and eventually we will take the long- wavelength, low-frequency limit w, fc — > 0. 

The form of the free fermion spectral function pp suggests that the above result can be written as a sum of the 
products p+ p- p+ p- and p- p+, where p± are the respective spectral functions for free fermion and 
antifermion [see eq. (|A2(l ]. The integrals over the dispersive variables po and So can be done trivially using the 
spectral functions for free fermions. For hard external fermion p ^ T and soft external photon k <^ T we can make 
the following kinematic approximations 

|p + k + q| - |p + q| ~ • k ^ fc , |p + k + q| + |p + q| ~ 2|p -I- q| ^ T . (3.27) 

After the analytic continuation of the external photon frequency ivk w + iO with to ^ k <^T, wc find the products 
P±P± lead to denominators of the form 
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whereas the products p±pz^ lead to denominators of the form 



1 



1 



2|p + q|Tc^ T 



(3.29) 



Hence in the fong- wavelength, low-frequency (or hydrodynamic) limit of the external photon fc, w — > 0, the products 
P±P± lead to singular denominators (pinch singularities) which furnish the leading contributions, whereas products 
of the form p±p^ lead to contributions that are suppressed by inverse powers of the temperature. We note that in 
our notation the pole-like singularities of the form given by eq. 1)3.28(1 must be understood with a principal value 
prescription as discussed above. Since these poles emerge in the hydrodynamic limit and play an important role in 
our discussion below, hereafter we will refer to them as hydrodynamic poles. 

Keeping only terms with the products p±p± that eventually lead to hydrodynamic poles, we find after some algebra 
that T^P'f^(^P -j- K, P) in the hydrodynamic limit can be written in the spectral representation as 



+ 00 



r^(so,P;it'fc,k) r^(so,P + k;wfe,k) 



So - «w„ 



So - ^tJJm - IVk 



where 



r''(so,p;fco,k) 7_(p) r+(so,p;fco,k) +7+(p) ^-(so, P; fco, k) 



(3.30) 



(3.31) 



The components F± are given by 



r±(so,p;ii^/c,k) 



q 



dqi 



it(p + q) 



(27r)3 

X <5(|p + q| - go - So) - 



p + q • k - zj/fc 

i-^(p + q) 
pTq • k + Wfe 

X ^(|p + q| - go + So) } [n.s(go) + 't-f(|p 



[Ktip, q) *PL{qO:q) + K^{P: q) Wi^O, q)] 
[K^iP,(l) *PLiqo,q) + Kf{p,q) *PT(go,g)] 

q|)] , (3.32) 



with L±{p) = (l,±p) being lightlike four-vectors. 

The soft-fermion contribution to the vertex has been studied in detail in Ref . . The conclusion of the detailed 
study of the vertex is that whereas the soft-fermion contribution is needed to fulfill the Ward identity, its contribution 
to the polarization is subleading Anticipating that in agreement with this conclusion that the soft-fermion 

contribution to the vertex will not contribute at leading logarithmic order, we will neglect this contribution for the 
moment and postpone a detailed discussion until we study the vertex corrections to the resummed two-loop photon 
polarization (see Sec. IIV 



E. Resummed one-loop Ward identity in the hydrodynamic limit 



In the program that we advocate here, namely the implementation of a resummation of the secular terms via the 
dynamical renormalization group in real time, we must confirm that the main ingredients in such a resummation fulfill 
the Ward identity, which guarantees that the result of the resummation will be gauge invariant. 

While eqs. (|3.22|1 - H3.23|1 and H3.24|l assert the fulfillment of the Ward identities, our focus below will be to extract 
the secular terms that are associated with hydrodynamic poles of the form (k • p — u!)~" in the limit k; w ^ 0. These 
poles are the manifestation of the pinch singularities and in real time they lead to secular terms in the perturbative 
expansion as discussed above. To lowest order with HTL resummed propagators these poles correspond to n = 1, 2. As 
it will become clear in the detailed discussion below, the dynamical renormalization group will provide a resummation 
of the leading order secular terms. Thus the building blocks of the resummation program are precisely these singular 
terms arising from the insertions of one loop self-energy and vertex with HTL resummed propagators. 

Thus we must confirm that the Ward identity between the resummed one-loop self-energy and vertex is fulfilled for 
these singular contributions. Once the Ward identity of these building blocks is confirmed, the dynamical renormal- 
ization group program that leads to a resummation of these leading order singular contribution is gauge invariant. 

As per the discussion above we will neglect the soft-fermion contribution to the vertex in the analysis that follows 
since its contribution to the polarization is subleading 0| (see below). The spectral representation of ^^P'^^ given by 
eqs. (|3.3U|) - ()3.32|) is particularly useful to establish the Ward identity between the resummed one-loop self-energy and 
vertex in the hydrodynamic limit. With the spectral representation for the self-energy eq. I|3.9|l in terms of ImE^'P 
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FIG. 3: Resummed two-loop photon polarization with the fermion self-energy correction. 



given by eq. (|3.11|l and the above spectral representation for the vertex F^'P'^, it is straightforward to check that the 
resummed one-loop Ward identity 



/ + C>C 7 
— ^ 
-oo 



ImE''P(so,P + k) ImE^P(so,p) 



So - lUJm - il^k 



So - it^m 



(3.33) 



is fulfilled in the hydrodynamic limit. The next step in the program is to insert the resummed one- loop self-energy 
and vertex in the transverse photon polarization, the fulfillment of the Ward identity between the self-energy and the 
vertex will ensure current conservation and gauge invariance. 



IV. RESUMMED TWO-LOOP TRANSVERSE PHOTON POLARIZATION 



At two-loop order the calculation of the photon self-energy involves two kinds of topologically different diagrams, 
namely the self-energy and vertex correction diagrams as depicted in Figs. |31 and ^ respectively. Similar diagrams 
have been studied for a hot QCD plasma within the context of photon or dilepton production (3Qi | and more recently in 
Ref. p9| but in a very different kinematic region. The color algebra, which introduces important differences, and the 
very different kinematic region of interest for dilepton production, which is not the relevant for transport phenomena, 
make these previous calculations not useful for the purpose of studying transport coefficients in the hot QED plasma. 



A. The self-energy contribution 



The two-loop diagrams for the fermion self-energy contribution to the transverse polarization are depicted in Fig. O 
In the imaginary-time formalism the self-energy correction to the polarization is given by 

nSE.«j(^) = e^TV Tr / [y E(P) S{P) 7^' S{P+K)+Y SiP) 7^ S{P+K) E{P+K) SiP+K)] , (4.1) 

where E(P) is one- loop resummed fermion self-energy in the spectral representation given by eq. H3.9|l and for later 
convenience we have written down explicitly the self-energy insertion for each fermion line. The Matsubara sum over 
the internal fermion frequency can be performed easily following the method illustrated in the previous section. After 
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FIG. 4: Resummed two-loop photon polarization with the fermion-photon vertex correction. 



using the identity eq. (|3.25|) we obtain 

.2 r r+oo 



^ -^li^L ""L ""L *{ 



X pF{to,p + k) 



dto<^ 7' Pf(po,p) IniI](go,p) Pf(so,p) 1^ 



+ 



So 



(qo - po)(go - so){to -qo- wk 



X pF{so,p + k) 



+ 7' PFito,p) 7^ PF(po,P + k) IniE((7o,P 
71^(^0) - npiso) 



k) 



(po - so)(po - go)(*o 



- Pa + i^k) (so - Po)(so - go)(to - so + ii^fe) 



(Q'o - Po){qo - so){tQ - go + i^k) 



(4.2) 



where ImI](goj?') is given by eq. (|3.1U|) and all the frequency denominators are understood with the principal value 
prescription. 

The imaginary part Imll^^'*^ (w, k) is obtained through the analytical continuation ivk — > w + zO and can be 
interpreted as cutting the photon self-energy diagrams in Fig. |3| in all possible manners. However, just as in the 
discussion of the vertex above, terms with products of the form p±p±p± give the leading contributions that feature 
hydrodynamic poles (or pinch denominators) in the hydrodynamic limit, while terms that involve mixed products 
of p± will be subleading by powers of 1/T. Therefore, only the products with the same particle or antiparticle 
components in all fermionic lines will yield the leading contribution. 

Focusing only on the leading contributions arising from terms where the fermion spectral functions are either all 
particle (p+) or all antiparticle (p-) [see eq. (|A2|) ]. we find 



ImnSE.y(^,fc) 



Tre^ PV 



(2^ 



dpo 



dqo 



+ 00 



dso 



X Tr [7' pf(po,p) Imi;((7o,_p) Pf(so,p) 1^ PF(io,P + k) 5{tQ - go 
-7* PF(io,p) 7^ PF(?3o,P + k) Imi;(qo,P + fc) pF(so,P + k) 5{to 



{qo-po)iqo - So) 
qo + (^)] 



(4.3) 



where the fermion spectral functions are either all particle or antiparticle contributions. Since the dispersive variable 
go is constrained by the delta functions S{to — qoT^), the leading terms that feature hydrodynamic poles are obtained 
from cutting the polarization diagrams in Fig. |31 through the fermion self-energy and hence the internal soft photon 
line. The other possible cuts do not lead to hydrodynamic poles, hence no secular terms arise from the other cuts. 
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FIG. 5; Diagrams corresponding to the resummed two-loop photon polarization in Fig. Qlb). 



Performing the integrals over pa, go, so, and to by using the delta functions, we find the leading terms in fc, — s- 
limit to be given by 



Imn^^'^Hto, fc) = 27r PV / 



(Pp dn p [p) 



f dp 



f - p^ 



Iml]+(|p + k| +ImS+(p + a>,|p + k|) 

(k-p-a;)2 



(4.4) 

where we have used the property Iml]_(— tjjp) — Iml]-|_(ti;,p). The above expression shows clearly that in fc, oj — > 
limit the leading contributions to Imll^^'*^ [uj, k) arise from terms that feature hydrodynamic poles of the form 



(k-pT^) 



(4.5) 



which will lead to linear secular terms in time. For k,Lu <ti p T we see that the imaginary parts of the fermion 
self-energy in eq. H4.4|l above are nearly on the mass shell of fermions or antifermions, the off-shellness being of order 
Iw — k • p|. Upon inserting Iml]+ given by eq. (|3.11l) into eq. (|4.4|l and noticing that for g <C T the second delta 
functions in eq. (|3.11|) cannot be satisfied near the fermion mass shell, we obtain 

ImnSE^^^(c.,fc) = 2vrV^PV J f fP J ^ J^^ dq, | [^^(go) + -i^dp + q|)] 



[K+{p, q) *pL{qo, q) + k:^ (p, q) *PT(go, q)] 
[Kt (p, q) *p+ (90, q) + (p, q) *p- (go, q)] 



1 



[I + nB{\p - (\\) - nF{qo)] 



2|p-q| 

(^(P ■ q - go - k ■ p + cj) + J(p ■ q - go + k ■ p - cj; 

(k.p-C.)2 



(4.6) 

+ (cj ^ —uo) . 



As will be shown below, the hydrodynamic poles in Imll^^'*^ (w, fc) correspond to secular terms in time in the function 
Q^^ [t) and signal the need for resummation in real time. 



B. Anomalous secular term from the self-energy 



At this stage we can begin establishing a connection with the program of resummation in real time by analyzing the 
leading secular terms arising in the perturbative expansion. As analyzed in detail in Sec. lIII Cl the leading contribution 
to the imaginary part of the self-energy near the fermion (or antifermion) mass shell is determined by the exchange 
of an ultrasoft transverse photon with momentum q <^ eT. Keeping only this leading contribution for the moment, 
we find that for fc <C p [see eq. (|3.16|) ] 



ImE+P(|p + k| -w,p) ImI]!P(p + cj, |p + k|) -aT In 



kp 



(4.7) 
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Introducing this leading estimate for the imaginary part of the fermion (antifermion) self-energy into eq. H4.4() . we 
obtain 



(Pp dn p {p) 
(27r)3 dp 



P'P^ In 



{uj -k-p) 



w — k • p 



(4.8) 



We now compute the function Q^^ (fc, t) by performing the Fourier transform in lo as per eq. (|2.10|1 . Using the properties 
of the principal value, we find 



1 - cos(yi) « (4.9) 



where Q^'^^''^^t) is given by eq. H3.3|l . With the help of the results established in Refs. [23,|23, we find the asymptotic 
long time limit to be given by 

gSE,<j (^) ^ giiTLaj ]^^2aTt\n{LJDt) + non secular terms] . (4.10) 

Therefore combining the one-loop (HTL) result with the leading contribution from ultrasoft transverse photon ex- 
change in the resummed self-energy correction to the two-loop polarization, we find 

g'^{t) = ^HTL.jj j-^ _ 2aTtln{ujDt) + non secular terms] . (4.11) 
Let us assume for a moment that the exponentiation of the leading secular term in the form 

9 

and explore the consequences of this result. [We used here eq. (|3.3|l ]. Such exponentiation would be justified consid- 
ering only the self-energy insertion on both fermionic lines in the loop and on the basis of the resummed form of the 
propagator eq. 1)3.17(1 l23,|23| since the function Q{k,t) behaves as the square of the propagator. Furthermore, as it 
has become customary in the literature, let us approximate 

^-2aTt\n{uJo t) _ ^~2Tt ^ (4.13) 

with the anomalous damping rate given by eq. 1(3. 18() . The time integral of G^-' (t) [see eq. (|2.30ll ] to obtain the 
conductivity can now be done in closed form and we find 

""TsF =YRi7;)+^"""s- 

This discussion highlights several important points: 

(i) In the limit /c, w — > the hydrodynamic poles of the form l/(w =F k • p)^, which are a consequence of the pinch 
terms |T3 |. result in secular perturbations for Q(k,t), i.e., terms that grow in time. The relation between pinch 
singularities and secular terms in perturbativc theory and their rcsummation via the dynamical renormalization 
group has been previously discussed in Refs. 22. 23^ ■ 

^From eq. (|4.9|l and the results obtained in Refs. H^ll^l we see that the hydrodynamic pole l/(a; — k-p)^ leads 
to a secular term linear in t, while the threshold singularity ln|cj — k • p| leads to the ln(w£i i) enhancement 
of the secular term. Thus, the extra logarithmic contribution to the secular term originates in the ultrasoft 
momentum region q <C eT of the exchanged transverse photon, which is only dynamically screened by Landau 
damping. Without this infrared, logarithmic enhancement the hydrodynamic pole would lead to a secular term 
linear in time. This observation will be important when we combine the contributions from the self-energy and 
vertex. 

(ii) // the imaginary part of the fermion self-energy on the mass shell were a finite constant F ^ T a Dyson 
resummation of self-energy insertions would lead to a Breit-Wigner form for the fermion and antifermion spectral 
functions near the mass shell, namely 

P^(Po^P) = l ^p^^p)2 + r^ - (4-15) 
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Using these spectral functions to compute the imaginary part of the polarization to one-loop order, keeping 
only the products p± p± and taking the limits fc, w — *■ 0, it is straightforward to find that the conductivity is 
independent of these limits and given by 

This simple exercise shows that the resummation of the secular terms in eq. (|4.12() by taking the coefficient of 
f to be a constant 2T gives the correct answer for the conductivity in the case where the damping rate is a 
finite constant and the spectral function near the quasiparticle poles is of the Breit-Wigner form, providing a 
reassuring check in a simpler case. 

(iii) In order to avoid confusion, we want to stress that the eqs. (|4.14|) - (|4.16|l are a result of the self-energy correction 
only. As will be discussed in detail below, the vertex correction cancels the ultrasoft photon exchange contri- 
bution from the self-energy and as a result the eqs. (|4.14|) - (|4.16|l do not apply to the hot QED plasma. The 
correct expression is given by eq. H5.89|l . The main point of the derivation of eqs. H4.14|l - (|4.16() is the following: 
(a) to illustrate what would be the result if only the self-energy correction is taken into account but without 
the vertex correction, (b) as we will show below that this approximation is equivalent to the relaxation time 
approximation in the Boltzmann equation. 

As it will be seen in detail later when we implement the dynamical renormalization group the resummed form of 
G(t) with only self-energy correction is much more subtle than the simple exponentiation assumed above in eq. 14.12|l . 
While the result obtained above is correct for a constant damping rate, the anomalous logarithmic time dependence 
will prevent the existence of a drift current at long times leading to a vanishing conductivity from only self-energy 
corrections from ultrasoft photon exchange. We will discuss these issues in more detail in Sec. lVBl where we introduce 
the dynamical renormalization group equation in the relaxation time approximation. 

While this discussion has highlighted these important points, keeping only the self-energy corrections is not con- 
sistent with the Ward identities and the vertex correction must be included. As advanced in Ref. ^25^1 and discussed 
explicitly below in Sec. IIVDI the contribution from the anomalous damping to the conductivity is in fact cancelled 
by the vertex correction. This cancellation, which will also be made manifest in the Boltzmann equation below is the 
reason that the conductivity is determined not by the quasiparticle relaxation time scale but by the transport time 
scale. 



C. The vertex contribution 



The vertex corrections to the photon polarization are displayed in Fig. ^ We focus first on the vertex correction 
with hard thermal loop resummed photon exchange, displayed in Fig.QIa). 

In the imaginary-time formulation the contribution to the polarization from this diagram is given by 



(K) = e'T ^ Tr 1 ^ 7' S{P) P (P + K,P)S{P + K), 



(4.17) 



where (P + K, P) is one- loop resummed fermion-photon vertex in the spectral representation given by eq. H3.30|l . 
The Matsubara sum over the internal fermion frequency is facilitated by the dispersive representation of the vertex 
r^{P + K, P) and of the fermion propagators. After using the summation formula eq. IjAlOjl and the identity eq. (|3.25|) 
repeatedly to combine factors, we obtain 



-e^Tr 



d^p 



+00 



dpo 



dSr 



(27r)3 

to - So - ivk to- po- ivk 

npito) - npipa) , npipo) - npiso) 



tn 



Po - il^k 



So - Po - ii^k 



dtf) Y Pf{po,p) 

^■^(s^,p;^^'fc,k) ^ 
So - Po 
r-^(so,P + k;wfc,k) 
So - to 



(4.18) 



/9F(to,P + k) 



where (so, P; ii^fe, k) is given by eqs. H3.31I) and H3.32|l and all the frequency denominators without ivk should be 
understood with the principal value prescription. 



21 



Upon the analytic continuation ivk lo + iO, the imaginary parts of the polarization arise from the following 
contributions: (i) The imaginary part of the denominator (delta functions) multiplies the real part of the function 
F-' ■ (ii) The real part of the denominator multiplies the imaginary part of the function F-' ■ We begin by analyzing 
the first contribution, namely those arising from the imaginary part of the denominators in the square brackets in 
eq. (HUl): 

(i) Terms proportional to 6{to — sq — uj) and 6{so — po ~ uj). These terms arise from cutting the polarization 
diagram in Fig. E^a) through the internal soft photon line because the dispersive variable Sq is constrained. For 
u <^ T the delta functions will have support only for the products p± p± . For such products of fermion spectral 
functions p±p± the denominator sq ~ po — tg — po — to and sq — to = Po ~ + ^ lead to the hydrodynamic pole 
of the form l/(tj=Fk-p) since po — and to — ±|p + k| for these products of fermion (or antifermion) spectral 
functions. Thus these discontinuities will lead to secular terms. 

(ii) Terms proportional to S{tQ — po — uj). These term arise from cutting the polarization diagram in Fig.^a) only 
through the internal fermion lines because the dispersive variable sq is not constrained. Again the delta function 
will have support only for the products p± p± in which case the contribution of this cut is of the form 



dp 



Rer^(so,P;w,k) ReF-'Xso,? + k; w,k) 



So TP 



soT |p + k| 



(4.19) 



where we have integrated over dispersive variables po and t^. Since the only hydrodynamic pole appears in 
the function F and is a single pole, these discontinuities will not lead to secular terms. Furthermore, upon 
relabelling p ^ p + k in the first term in the momentum integral, we find that these contributions cancel each 
other. 

We now analyze the second contribution, namely those arisiiig_from of the imaginaryjDart of the function F-' . The 
imaginary part of the function F-' is obtained by replacing p + q • k =F fco — > ±i7r(5(p + q • k ^ w). The terms that 
can give rise to hydrodynamic poles of the form l/(k • p =p w) are the second and third terms in eq. (|4.18|l when the 
product of fermion spectral functions is either p+ or p_ p_. Combining these two terms we obtain contributions 
of the form 



dnpip) k • p 
dp =F k • p 



ImF-' (so, P; w, k) ImF-' (so, p + k; fco, k) 



So- P 



So 



|P + k| 



(4.20) 



There is a single hydrodynamic denominator in this expression since the hydrodynamic denominator in F-' became a 
delta function. Hence there is no secular term arising from this discontinuity. Furthermore, upon relabelling p — > p+k 
with fc <C p in the first term in the momentum integral, we find that these contributions vanish and hence they do 
not contribute in the limit fc ^ 0. 

After this analysis, we conclude that the leading contributions from the vertex correction in Fig. 0Ja) that feature 
hydrodynamic poles leading to secular terms and therefore contribute to the conductivity are given by 



npito) - "-f(so) 
So - Po 



ne^ Tr 



d^p 
(2^ 



dpo 



dsQ 



dto 7' PF (po , P) (^0 , P + k) 



(4.21) 



Rer^(so, p;ll>, k) S{to - sq - lo) + ^^^^'^^ npiso) j^g/rJ^^^^ p + k; k) S{so - po ~ u) 

So - to 



where ReF-' (sq, p; w, k) can be extracted from eqs. H3.31|l and (|3.32l) through the analytical continuation iiy^ —^ui + iO 
and all frequency denominators [including those in ReF-'(so, p; w, k)] should be understood as the principal values. 
Furthermore, the fermionic spectral functions in the above expression are both p+ or both p_ since only these products 
lead to the hydrodynamic poles. 

Using the delta functions to perform the integrals over the dispersive variables sg, po-, and to, we find the leading 
terms in fc, w — > limit to be given by 



lmUy^{uj,k) = -2^e2a;Pv|^^^|j-|— [Reri(|p + k|-c.,p;c.,k) 
+ ReF^ip + p + k; tj, k)] + {lj -u 



(4.22) 
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where we have used the property ReF_(— so,p;ti;,k) = ReF+(so, p; — w, k). Upon inserting ReF+ obtained from 
eqs. and H3.32|l and keeping in ReF+ only the delta functions that can be satisfied for cj, fc ^ 0, we find 

lmny\co,k) = -2.e4c.PV J J ^ /^J d<Zo [^^(go) + n^(|p + q|)] 



(27r)3 dp J (27r)3 

X f f+~q [K+{p, q) *pL{qa,q) + K+{p, q) *pT{qa, q)] 
'S{p • q - 90 - k • p + cj) + 5(p • q - 90 + k • p - cj) 



+ ^ -uj) 



(4.23) 



(k • p - cj)(k • p + q - cj) 

The contribution from the diagram in Fig. 01b) to the resummed two- loop photon polarization corresponds to those 
from the two diagrams displayed in Fig. |31 which can be now calculated handily by writing each one-loop vertex in 
terms of the spectral representation similar to that given by eq. (|3.30() but with the free photon spectral functions 
replacing those in eq. (|3.32|) . In order to simplify the computation we will set the external photon momentum fc = 0. 
The sum over the Matsubara frequency in the soft fermion loop can be easily done. The resulting expression is rather 
lengthy but the imaginary part can be computed straightforwardly following the analysis presented above. In the 
limit ^ 0, we find 



Imn^'*^(w,fc^O) 



d^p 
(2^ 



1 



1 



(5(mo - 



- LO 
So) 



to - pa+u 
dnpiso) 



dpn 

-OO 

dnpiso) 
dsQ 
1 



dSr 



+ OC 



dta ReF!,(so,P) *P+{Po,p) *P+{to,p) 



[SiPo - So) + S{to - So)] + 2 
1 1 



duo ReF+(wo, p) 



(io - so)(Po - So) dso uj{uo-sa)'^ \to - po - h - po + 
[o[Po - So) + 6{po - uo)\ H [d[to - so) + d(io - uo)J 



dpo 



dto 



(4.24) 



where F±(so,p) = a; F±(so, p; t^^, 0) and F±(so, p; w, 0) is given by eq. (|3.32(l but with the replacement *pT PB 
and *pL 0. In obtaining the above expression, we have used the property F-L(— so,p) = — F^(so,p), *p-{—pQ,p) — 
*p+{po,p) and safely set w = inside the arguments of the delta functions and in the denominators which do not 
vanish in this limit. 

Let us first focus on terms proportional to I /{to — po T t^) in cq. (|4.24(l . Whereas there are pole-pole, pole-cut 
and cut-cut contributions in the product *p+{to,p) *P+{potP)^ only the pole-pole contribution will lead to pinch 
denominators. However, upon inspecting the support of the delta functions in F [see eq. (|3.32ll and recall the above- 
mentioned replacement] , it becomes clear that the pole-pole contribution vanishes because the quasiparticle pole aj+ (p) 
lies above the light cone but the delta functions in F only have support below the light cone. Thus only the terms 
proportional to feature pinch denominator and hence lead to linear secular terms in Q{t). However, will argue 
below that these linear secular terms are subleading in the leading logarithmic approximation. Therefore, we conclude 
that while the contribution from the diagram in Fig.01[b) is necessary to fulfill the Ward identity, it actually does not 
contribute to leading logarithmic order in agreement with the conclusion of Ref. 49]. 



D. The cancellation between vertex and self-energy leading contributions and the extraction of leading 

secular terms 



Gathering the above results for the self-energy and vertex contributions [see eqs. (|4.7|l and (|4.23() ]. we find that the 
resummed two-loop photon polarization Imn*-'(cL;, k) can be written in a compact form as 

Imn2-i°°P''J (w, /c) = ImH^'^''^' {u,k) + ImH^^^^' {uo, k) 

d^p dnpip) f d^q '■+°° 



-27re^cj PV 



(27r)3 dp J (27r 



dqo 



[r\Lo,k,p,q)-r\Lo,k,p,o)] 



["-^(go) +"-F(lP + ql)] [i^+(p,q) *pLiqo,q) + K+{p,q) *pT{qo,q) 



2lp-ql 



X [1 + nB{\p - q\) - nF{qo)][Kt {p^q) *p+{qo,q) + (p,q) *p_(9o,9)] 
X [5{p • q - go - k • p + cj) + (5(p • q - go + k • p - tj)] + (lj ^ -uj) \ , 



(4.25) 
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where 

^^■'■(a;,k,p,q) = ^ P' P + g^ ^ (4 26) 

(k • p - tj) (k • p + q - tj) 

This expression makes clear that the region of ultrasoft photon momentum q <^ eT p ^ T which leads to the 
anomalous damping eq. (|3.18|) is cancelled between the contributions of the self-energy and those of the vertex to 
the polarization. This is the cancellation that was found in Ref. psf . which our analysis makes manifest in the 
perturbative framework. Whereas the region of ultrasoft photon momentum is suppressed in the electric conductivity 
as a result of the explicit cancellation between the self-energy and vertex corrections, there will be contributions to 
the conductivity from the region of exchanged momentum eT q T ^ which within our real-time framework will 
arise as linear secular terms in the perturbative expansion. 

The analysis of the secular terms in the polarization studied in the case of the self-energy correction in Sec. IIVBI 
revealed that the hydrodynamic poles (or pinch denominators) of the form In jcj — k • p|/(aj — k • p)^ give rise to a 
secular term of the form tint in Q{t). The logarithmic enhancement of such secular term is due to the logarithmic, 
infrared divergence in the self-energy that arises from the region of ultrasoft photon exchange q <^ eT. However, since 
the soft-photon contribution in eq. (|4.25() involves the difference of the function which vanishes for q = 0, for q <^ p 
we can expand the integrand in powers of q/p. Rotational invariance dictates that only even powers will survive the 
angular integration. Thus, the lowest order term in the g-integral will have an extra power of q^ in the numerator 
which will render finite any potential logarithmic, infrared divergence that is responsible for the anomalous damping. 
This extra power of q^ in the numerator is a manifestation of the fact that the electrical conductivity is determined 
by the transport time scale @. 

We can now proceed to calculate the corresponding Q{k, t) function in the limit fc ^ 0. Since Imn^"'°°P'*^ (w, k) is 
free of logarithmic, infrared singularity responsible for the anomalous damping, the Fourier transform in eq. (|2.10f) 
can be performed easily by using contour integration. In particular, we need 

PY d:^- c-''^' ^(P-q-goTk-p±cj) _ iTrS{p-q-qo) ^^-^k-pt ^-»k.^tj 
J -00 (k • p - w)(k • p+~q - w) k-(p'+q-p) 

-TTte-^'^-P* J(p-q-qo) , (4.27) 

where in the second approximation we have extracted the leading secular term in the small fc, q limit. We find to 
two-loop order with hard thermal loop propagators and in fc ^ limit 

f-'°°^^'^t) = -8neH [ / dq, Up^fP-ffTq') K(go) + n^^dp + q|)] 

—00 L 



(27r)3 dp J (27r)3 



[K+{p,c\) *pL{qo,q) + K+{p,q) *PT{qo,q)] + —r—AK+{p,q) *p+{qQ, q) + {p,q) *p-{qo,q)] 

^ IP Ql 

X [l + nB(|p-q|)-nF(go)]| %o - P • q) • (4.28) 

For soft momentum q ^ p we can expand the integrand in powers of q/p. The cancellation between the self-energy and 
vertex for ultrasoft photon momentum transfer q thus entails an extra power of q^ in the soft-photon contribution 
to the integrand in eq. (|4.28|l . which as remarked above is the origin of the transport time scale. 

To compute a from eq. H2.3()|l we ultimately need Q{t) = We can obtain a rehable estimate for the 

leading logarithmic contribution by utilizing the following approximations. Since *pl,t{(1q, q) are odd functions in go, 
we approximate 

nB{qo) + nF{\p + q\) ^T/qo , 1 + ubHp - q\) - npiqo) ^ 1/2 + nsip) . (4.29) 

Furthermore, following Baym et al. f^, we neglect the HTL contributions in the denominators of the respective HTL 
spectral functions 



1 ImngTL(g,go) e^T^ qo „ , 2 2^ 

p^(io,<i) - -- — — = ^7®^'' "''^ 



1 ImnHTL(g,go) e^T^ qo 1 2 2 

TT [q^ - q^Y 12g4 g l _ % 



^ iq,^q)2 - T6^ ^('^ --^o) ■ (4.30) 
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FIG. 6: Three-loop diagrams for the photon polarization in the leading logarithmic approximation. 

Performing the trivial qo integral in eq. 14.28|l and keeping the leading terms in q/p in the integrand, we find that 
the leading logarithms arise from the momentum integral of the form J^"'"' dq/q = ln((7max/'Zmin), with the upper and 
lower momentum cutoffs gmax ~ T and Qmin ~ eT, respectively, which limit the validity of the approximations made. 
We thus see that the leading logarithmic terms emerge from the region of exchanged momentum qmin ^5 9 ^ 9max and 
appears in the form ln((7inax/Qmin) ~ ln(l/e). In particular, this entails that in order to extract the leading logarithms 
we can simply use the perturbative spectral functions. 
Finally, upon collecting the results we find 

g'-'°°^t) = -^^^ t;"^^ T In t , (4.31) 

where [see eqs. and (|2.29ll ] 

gKTL^_^^ (4.32) 
lo 

Therefore, up to three-loop order we obtain that 

(4.33) 

This expression clearly displays the linear secular term arising from denominators of the form l/(w — k • p)^ without 
the logarithm of time featured by the anomalous damping precisely because of the cancellation between self-energy 
and vertex for ultrasoft photon momentum transfer. It also highlights that the transport time scale is given by 
tti ~ [e'*Tln(l/e)]~^, which is much longer than the fermionic quasiparticle relaxation time scale tqp ~ [e^Tln(l/e)]~^. 
The origin of the logarithms in these expressions is also different: in the transport time scale the logarithm originates 
from momenta eT ^ q T, while in the quasiparticle relaxation time scale it is momenta e^T ^ 9 ^ eT that lead to 
the logarithm. Our analysis reveals that this difference is a consequence of the cancellation between self-energy and 
vertex corrections for ultrasoft photon momentum transfer. 



g{t) = g 
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FIG. 7: Cuts in the three-loop diagrams for the photon polarization in the leading logarithmic approximation that lead to the 
secular term. The dashed lines determine the cuts. 



Origin of the leading logarithms 

An important aspect gleaned from the above result is that the leading logarithmic correction can be extracted by 
using the perturbative spectral functions given by eq. H4.30|l and by cutting off the integral over the transfer momentum 
q with upper and lower cutoffs (/max ^ T and Qmin ^ eT, respectively. This simplification corresponds to computing 
the three-loop diagrams depicted in Fig. with cutoff momentum integrals. The imaginary parts of the polarization 
that lead to the secular term and leading logarithmic contribution correspond to the cuts displayed in Fig.[7| The cut- 
diagrams with the photon self-energy correspond to M0ller scattering — > e^, since the singular denominators 
that lead to secular terms arise from the products and the photon spectral function below the light cone 

in the HTL approximation arises also from these type of products in the photon polarization bubble. The internal 
longitudinal (transverse) photon propagators with momentum q contributes the factor (g^)^^ [{qq — 9^)"^] that leads 
to the logarithms. The cut-diagram with the fermion self-energy corresponds to Compton scattering j ^ j as 
well as pair annihilation and creation processes e"*" e~ ■'-^ 7 + 7, where the electrons, positrons and photons are all 
hard. The internal fermion propagators with momentum q give rise the factor (go T that leads to the logarithms. 

The diagrams that yield the leading logarithms in the coupling are formally corrections of 0{e^) to the photon 
polarization. The logarithms arise from the free intermediate state propagators with momentum restricted by the 
upper and lower cutoffs gmax ^ T and qmin ^ e T, respectively. Thus, only the polarization diagrams depicted in Fig.d 
with the simplified spectral functions eq. H4.30|l and the upper and lower momentum cutoffs qmax ^ T, qmin ^ e T 
(or, alternatively, only the cut diagrams shown in Fig. with restricted momentum transfer eT ^ q '^T) contribute 
linear secular terms to leadinglogarithmic order. Such simplification in the leading logarithmic approximation has 
already been noticed in Ref. j5(. In turn, this implies that both Debye screening for the longitudinal photon and 
dynamical screening by Landau damping for the transverse photon propagators are irrelevant to extract the leading 
logarithmic contribution and we can simply extract the leading secular terms to leading logarithmic accuracy by 
computing the three-loop corrections to the polarization with free vacuum photon propagators and upper and lower 
momentum cutoffs of order T and eT, respectively. This observation will be important below where we establish the 
correspondence with the Boltzmann equation. 

This analysis also allows to identify the diagrams in the perturbative expansion that yield the leading logarithmic 
contribution and the linear secular terms. To leading logarithmic order there are additional corrections to the polar- 
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FIG. 8: Three-loop photon polarization diagrams of 0{e^) that feature a linear secular term but no logarithms in the coupling. 

ization formally of 0(e®), as depicted in Fig. |S1 They correspond to a vertex correction to the fermion self-energy 
and a crossed ladder diagram. The calculation of these diagrams is very complicated and beyond the scope of this 
article. However, from the discussion above it is clear that these diagrams will feature a linear secular term but no 
logarithms of the coupling. 

This above discussion highlights the origin of the leading logarithmic contributions and we can now argue firmly 
that the contribution from the diagram in Fig.^b) (or, equivalently, those depicted in Fig.jSJ is subleading at leading 
logarithmic order. As argued after eq. 14.24(1 . it is the term proportional to l/oj^ that can lead to secular terms. To 
leading order we can replace the HTL fermion spectral function *p+{po,p) in eq. 14.24|l by that given in eq. H4.30|l . 
As a result, the contributions from the diagrams in Fig. El are of 0(e®) and ©(e^^), respectively. Therefore, the 
soft-fermion contribution to the vertex is subleading in agreement with the conclusion of Ref. jl^ . 



V. DYNAMICAL RENORMALIZATION GROUP AND BOLTZMANN EQUATION 

In the previous sections we have extracted the corrections to the photon polarization that feature linear secular 
terms in time to leading logarithmic approximation. The next task is to sum the secular terms in the perturba- 
tive expansion by using the dynamical renormalization group approach introduced in Ref. 12,'^ . The dynamical 
renormalization grou p appl ies to the equation of motion of expectation values, either of the Heisenberg fields or of 
composite operators [23, l23| which lead to an initial value problem. Thus the first step is to obtain the corresponding 
equations of motion for Q{t) or an alternative quantity from which it can be derived. In order to understand what 
are the equations of motion (and corresponding initial value problem) for Q{t) (or a related quantity), we begin with 
equation eq. H2.28|l which basically identifies G {t) with the time derivative of the expectation value of the current in 
the zero-momentum limit. Thus we now study the expectation value of the current to identify the proper equation of 
motion in order to implement the DRG resummation. 

We begin the analysis by writing 

*(x, t)=Y^j ^ [6p,A(<) «pa e-^^^'-P-) + dl^{t) «p,, e^(P*-P-)] , (5.1) 

where u and v are the free massless Dirac spinors. The time evolution of the operators b-p^\{t) and dp^x{t) is determined 
by the Heisenberg equations of motion. Here we have factored out the rapidly varying phases e^****, therefore for 
momentum p ~ T ^ m these operators vary slowly only due to the interaction. An important aspect of the gauge 
invariant formulation introduced earlier is that the operators b and d are manifestly gauge invariant, since the field 
operator ^E" has been constructed to be manifestly gauge invariant. 

Using the properties of the free massless Dirac spinors u and v, we find that the zero momentum limit of the 
expectation value of the electromagnetic current is given by 



(5.2) 
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with 



m ^ eY^J^f [{bl,{t) fep,,(t)) - (d_p,.(t) dl^Jt))] , (5.3) 

.m = -^eJ2 I 7S3 [(^LaW dlx'it)) (P ^^Y ^-p,A' e^'^*] + c.c. , (5.4) 

where are the spin matrices. The term 3c{t) is the convection current and 3s{t) is the spin current. The convection 
current term is time independent in the free theory and is therefore varying slowly in the interacting theory, the time 
variation is solely due to the interactions. On the other hand, the spin current term mixes particles and antiparticles 
and oscillates very fast (even in free field theory), on time scales ~ 1/T. This term is responsible for the 
phenomenon of zitterbewegung associated with the mixture of positive and negative energy solutions in relativistic 
fermionic wave packets. Thus under the assumption that processes associated with photon exchange will not mix 
particles and antiparticles, namely that the exchanged photon momenta are <ti T, we can neglect the contribution 
from the spin current term since it oscillates very fast and will average out on time scales ^ 1/T. 
Introducing the spatial Fourier transforms of the fcrmion fields 

vl/(x,t) = / *(P,i) e'P'' , ^H^,t) = / *^(P'^) ' (5.5) 

we find the following expressions for the expectation values that enter in the convection current 

E(^p,aW bpAt)) = (*(P,^) 7-(p) ^iP,t)) , 5](d-p,A(0 dLp,;,(t)) - (*(p,i) 7+(p) *(P,t)) • (5.6) 

A a 

These are identified as the contributions from positive and negative energy components to the convection current. 
Under the assumption that all the processes involved do not lead to a mixture of positive and negative energy 
components, the two contributions will not mix in any order of perturbation theory. This is manifest in the perturbative 
calculation of the previous section in that only products of spectral functions for particles or antiparticles enter in the 
secular terms. Thus the slowly varying part of the current at zero momentum is determined by the convection term 
and given by the expression 



' J^P' [(*(p,t)7-(p)*(p,i))-(*(p,i)7+(p)*(p,t))] =e J ^ f {^{p,t) j-p ^{p,t)) . (5.7) 
It proves convenient to introduce the spin-averaged expectation values of the number of particles and antiparticles as 

/(p, t) = l 5](6^,,(t) 6p,A(t)) , 7(P, t) = ^ E<<,A W dpAt)) , (5.8) 

A A 

which as emphasized above are gauge invariant. In terms of these, the convection part of the current is given by 

^2ej^p^ [/(P,t) +7(-p,t)] , (5.9) 

which is the same form as that in kinetic theory, hence we call this the kinetic current. This analysis clearly highlights 
two important aspects: (i) The kinetic current is identified with the expectation value of the convection current in 
Dirac theory, this is the component that varies slowly on the time scale l/T, the contribution from the spin current 
averages out on time scales larger than 1 /T provided there are no processes with transferred momenta ~ 2T that can 
mix particles with antiparticles. (ii) The distribution functions / and / are identified as the expectation values of the 
number operators for particles and antiparticles, respectively. Our gauge invariant formulation guarantees that these 
are indeed gauge invariant physical observables. 

^From now on we will refer to the current simply as the convection part and the discussion above clearly indicates 
that this is the kinetic current. Thus we can establish the following correspondence between the quantum field theory 
and the kinetic approach [see eq. H5.2|l ]: 

.P(t) = -ie j ^fTT[-fpS<{p;t,t)\ , 
/(p,i) = -^Tr[7-(p)5<(p;t,t)] , 

7(-p,t) = l + ^Tr[7+(p)^<(p;t,i)] . (5.10) 
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The correspondence given by eq. (|5.10() makes manifest that the kinetic current is determined by the equal-time hmit 
of the spatial Fourier transform of the full fcrmion propagator. The full fermion propagator is given by a Dyson sum 

5 = 5o + 5o S 50 + So S 5o SS-o + • • • . (5.11) 

If an external electric field S = —A-t is present, then So is the free propagator in the presence of the background 
gauge field At but without corrections from interaction due to the fluctuating gauge fields, which are contained in 
the self-energy E. This observation for the kinetic current makes explicit that the collision term that enters the 
Boltzmann equation is determined by the fermion self-energy (see below). 

Furthermore, at this stage we can make contact with the function Q^^ (t) defined by the linear response relation 
eq. (|2.28|l . For this purpose we define to linear order in the external electric field 

fip,t)=nFip)+£-pAip,t) , J{~p,t)^nF{p) + £-pA{p,t). (5.12) 

We therefore identify, 

g^'it) = -2e J ^ f p' [a(p, t) + t{p, t)] , (5.13) 

where the dot stands for derivative with respect to time. ea. H2.30|l then leads to the following expression for the 
conductivity 

2e f cPv — 
CT = — J [Aip,t ^ cx,) + A{p,t = ^)] , (5.14) 

where we have used the initial condition A(p, t = 0) — A{p, t = 0) = 0. 

We now derive the Boltzmann equation for the distribution functions / and / in the presence of the background 
gauge field by implementing the dynamical renormalization group. We first calculate the time derivative of these 
distribution functions to linear order in the background gauge field as in Sec. IIIBI and extract the positive and 
negative energy components of the coefficient of p^ from the HTL photon polarization. As discussed in Sec. IIII iO the 
polarization for small external frequency and momentum only receives contributions from the products of fermion 
spectral functions of the form and p_ p_ , corresponding to the positive and energy fermion mass shells 

respectively. This explicitly shows that there is no mixing between positive and negative energy components in the 
small transferred momentum limit, thus justifying the neglect of the spin current contribution. 

We assume that at a given time to the initial density matrix is diagonal in the number basis and that the number 
operators for fermion and antifermion have respective expectation values /(p,io) and /(p,to) as given by eq. H5.8|l . 
This choice of the initial density matrix entails that the fermion propagators have the free field form as given in 
Appendix IbI but in terms of the nonequilibrium distributions /(p, to) and /(p, to) at time to (see also Ref. |2^). From 
the results of Sec. lIII Al wc find to lowest order in the background field and leaving the calculation of the collision terms 
(arising from the interaction with the photon fields, namely the self-energy corrections) for the moment, that the time 
derivative of the distribution functions are determined by the HTL photon polarization. Using the nonequilibrium 
fermion propagators given in Appendix ^ we find 

d d - - 

Q^fiP^t) — —e£ •Vp/(p, to) -I- collision terms , — /(— p,t) = —e£ •Vp/(— p,to) + collision terms . (5.15) 

The next step is to obtain the contribution to the time derivative of the distribution functions from the interaction, 
namely the collision term. From the definition of the distribution functions eq. I|5.8f) . it follows that their time 
derivatives are completely determined by the Heisenberg equations for the fermion fields, namely, 

—■^{p, t) = -ia-p ^'(p, t)+ie j [a ■ AT(q, t) - A°(q, t)] *(p - q, t) . (5.16) 

Taking the time derivative of the distribution functions and using eq. H5.16|l . we find that the collision term for /(p, t) 
is given by 



9/(p,t) 



dt 



coll 



le 
"2 



/ [('I'- (p + q, 0^°'- (q, 07- (P)*+ (P, t)) 



-(*-(p + q,t)7.A^(q,t)7V(p)*+(P,i))] + cc. , (5.17) 
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FIG. 9: Diagrams corresponding to the collision term df{p,t)/dt\coii up to one-loop order: (a) the soft-photon contribution 
and (b) the soft-fermion contribution. The solid line denotes the free fermion propagator in the presence of the background 
gauge field. The heavy cross in the fermion line represents the insertion of 7^7- at the time t. 

where the superscripts '±' refer to fields in the forward and backward branches in closed-time-path formaUsm, re- 
spectively (see |23| for details) . The collision term for f{—p,t) has a similar expression but with 7-(p) — > — 7+(p). 
The expectation values in these expressions are in an initial density matrix and can be systematically computed in 
perturbation theory. In order to compute the time derivative of the distribution function in eq. H5.17|l we now follow 
the procedure outhned in Sec. IV B in Ref. [23l |. 

Up to one-loop order, as depicted in Fig. |51 the collision term receives contributions from soft internal photon and 
soft internal fermion 



Following the steps detailed in Ref. |23.l23l| with the propagators given in Appendix^ we find that the collision term 
can be written in an illuminating form which displays the familiar gain minus loss form as 



9/(p,<) 



to), /(., to), t^to]= to), /(., io), t - to] + C^'ifi; to), /(., to), t - to] . (5.18) 



dt 



coll 



C[/(.,to),/(.,to);t-to] = 2 



— 00 



dpo 



sin[(p-po)(t-to)] 
Tr{p - po) 



{[l-/(p,to)] S>(po,p)-/(p,to)E<(po,p)} , (5.19) 



where 



^>{Po,p) = ^>{pq,p) + S>(po,p) , ^<{Po,p) ^ ^<{Po,p) + S<(po,p) , 



(5.20) 
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with the following self-energy components 

"+00 



S>Pbo,p) - j J dqo{[K+{p,ci) *pLiqo,q) + K+ip,<i) *pTiqa,q)] 



/(p + q,to) nsiqa) S{po - go - |p + q|) + [K^ (p,q) *pL{qo,q) + (p, q) *PTi<la,q)] 
[1 -7(-p- q,io)] nB{qo) S{pa - qo + |p + q|)} , 
d^q '■+°° 



^<iPa,p) = e^TT J -^-^ J dqo {[K^{p,q) *pL{qQ,q) + K}{p,q) *pT{qo,q)] 



(27r)3 

[1 - /(p + q,to)][l + n-i3('Zo)] S{pa - go - |p + q|) + [-ftrr(p,q) *pLiqo, q) + {p,(i) *PTiqo,q)] 
x7(-p-q,to) [1 + nsiqa)] Sipo - qo + \p + q\)} , (5.21) 



X 



and 



d^o 1 



J (2^2|p~q| y {[-^^CP.q) *P+iqo,q)+K^ (p,q) VC^o,?)] 

n-F(go) J^sdp- q|) S{po -qo - |P- q|) + [/^f (p,q) *p+{qo,q) + K}ip,q) *p^{qo,q)] 
[1 - n F iqo)][l + riB i\p - q\)] 5{po + go + |p - q|)} , 
f d^q 1 f^°° 

J (2^ 2|p- q| J [^'^ ^) + ^3 (p. q) *p-{qo, q)] 

X [l-nF(go)][l + nB(|p-q|)] (5(po - 90 - |p - q|) 

+ [-ft-VlP^q) *P+(9o,9) + ^3^(p,q)*P-('?o,'?)] nF(go) n-B(|p-q|) (^bo + go + lp-ql)} . (5.22) 
For antiparticles we find that 

d7ip,t) 



, ,J3„ 1 r+00 

^%{po,p) = ne' 



C[/(.,to),/(.,to);t-to] , (5.23) 



dt coll 
where 

C[f{P. to), 7(-q, to); t - to] = C^[7(-P, to), /(q, to); t - io] , (5.24) 

which entails that /(p, t) and /(— p, i) satisfy the same evolution equations. Since the evolution equation for antipar- 
ticles (with momentum — p) is identical to that for particles (with momentum p), we now focus solely on the particle 
distribution function, the steps can be repeated straightforwardly for the distribution of antiparticles. 

A. DRG resummation of the secular terms 

Combining eq. H5.15I) with the collision term obtained above we find 



dt 



-/(p, t) = -e £ • Vp/(p, to) + C[/(., to), /(., to); t - to]. (5.25) 



Notice that the right-hand side depends only on time through the sine functions in the collision term, and the 
distribution functions /(p,to) and /(— p,to) are fixed at the initial time [see eq. 15.19|l ]. Thus the evolution equation 
can be integrated to yield 

/(p, t) = /(p, to) - e 5 • Vp/(p, to) (t - to) + f C[/(., to),7(., to); t' - to] dt' . (5.26) 

Jto 

The last term gives upon integration on time, 

f C[f(,,to)J{.,toy,t'-to]dt' ^ r°°rfpo l'^"^[(P-Po)(^-^o)] /C(/(.,to),7(-,to);Po,p), (5.27) 



with 

/C(/(.,to),7(.,to);Po,p) = 2{[l-/(p,to)] E>(po,p)-/(p,to)S<(po,p)} . (5.28) 
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In general the time integral of the collision term will feature secular terms, whose time dependence is determined by 
the behavior of the term that multiplies the cosine near the resonance p = pq. For example, if the coefficient of the 
cosine term IC{f{.,to),f{., to);po,p)is finite and does not feature any singularity as p — po ^ then we can use the 
long-time asymptotic expression '2^ 

l-COs[(p-po)(t-to)] /, , W/ N /r 0Q\ 

"2 > (t - to) d(p - po) (5.29) 

to extract the secular term, which in this case is linear in time as is the case in the usual Fermi's Golden rule. 
However, if /C(/(., io): fi-i ^o);PojP) either vanishes or features singularities near the resonance p = Po then one must 
be more careful in extracting the secular term. Such is the case with the anomalous logarithms that will emerge in 
the relaxation time approximation (see below). Whatever the form of these secular terms, their main feature is that 
they grow in time and invalidate the perturbative expansion, therefore they must be resummed. This resummation 
scheme in real time is the goal of the dynamical renormalization group program p^ . l2^ . 

Just as the usual renormalization group absorbs the divergences into a renormalization of couplings at an arbitrary 
momentum scale, the dynamical renormalization group absorbs the secular terms into a finite renormalization of the 
distribution function at an arbitrary time scale r. Let us define 

/(p,io) = /(p,r)Z(p,r,io), (5.30) 

where the renormalization constant 2 will be chosen in an expansion in £ and the coupling that enters in the collision 
term to cancel the secular terms in eq. (|5.26|l at the time t = t. Choosing to lowest order 

e£.Vp/(p,T),^ ^ , r CSi;r)J{.,Ty,t' -to] , ,^,^2 ^2 



z(p,T,to) = i+ ^^ ' {r-to)- 'y,'" — ^d^'+o(£^c^£:a), (5.31) 

where Cg denotes the part of the collision term that leads to secular terms. Then to lowest order in the external 
background field and in the strength of the collision term we find 

/(P,<) = /(p,T)-e5-Vp/(p,r)(t-T)-^ a[/(.,T),7(.,T);t'-to] dt'+0{£^,Cl£ C.)+non secular terms , (5.32) 

Choosing the arbitrary scale t near t, the perturbative expansion is now valid in a neighborhood of ^ = r by 
renormalizing the distribution function. 

However, the distribution function f{p,t) cannot depend on the arbitrary time scale r, namely df{p,t)/dT = 
leads to the dynamical renormalization group equation to this order (after setting t — t) 

^^f{p,t) + eS-Vpf{p,t) = Cs[f{.,t)J{.,t);t] . (5.33) 

This dynamical renormalization group equation is recognized as the Boltzmann equation, where the collision term is 
the part of the original one that leads to secular terms. (See Ref. _22] for the scalar theory). Thus in order to write down 
the Boltzmann equation we must first work out Cs- As we shall see below in sec. V C, when /C(/(., to), f{-,to)]Po,p) 
is regular at po = p, we simply find 

C«[/(.,t),7(.,i);i]=/C(/(.,t),7(.,i);p,p) . 



B. Relaxation time approximation 



In the relaxation time approximation only the particles with momentum p and the antiparticles with momentum 
— p are slightly distorted out of equilibrium, while particles and antiparticles with all other momenta are assumed to 
be in equilibrium. Hence, we write the particle distribution as 

f{pM)^nF{p) + 6f{p,to) /(p + q,to)=nF(|p + q|) for q^O, (5.34) 

and similarly for antiparticle distributions. With the above distributions the collision term given by eq. (|5.19|l simplifies 
to 

C[f{.M)J{.M)\t-to]^-2 5f{pM) r^dpo ^+{po,p) ^'""^(P - Po)(t - to)] ^ ^53^^ 

J -00 T^KP-Po) 
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where we have used the fact that when evaluated with the equihbrium distribution (npip)), }(^(nF{-),np(.)]pQ,p) 
vanishes by virtue of detailed balance and I]>(po,p) + ^<iPojP) — S+(po:P) with E+(po,p) given by eq. H3.11I) . 

As discussed in Sec. lIII Cl the self-energy E_|_(po,p) near the fermion mass shell is dominated b y th e ultrasoft photon 
exchange that leads to anomalous damping. Using eq. H3.16|l and the results quoted in Refs. "2^, 23] , the time integral 
in eq. H5.26|l to linear order in the external electric field and the departure of equilibrium (considered to be of the 
same order as the external electric field) is given by 

Sf{p, t) = <5/(p, h)-e£ ■VpnF{p){t - to) - 2aT {t - to) In [{t - to)LUD] Sfip, + subleading . (5.36) 

This expression clearly displays the anomalous logarithms as secular terms found previously in the perturbative 
expression eq. (14.10(1 . 

Following the steps described above by introducing the renormalization constant Z(p,T, tg): 

5/(p,io)=2(p,r,to)J/(p,r) , (5.37) 

with 

Z(p, T,to)^l + e£ ■VpnF{p){T - to) + 2aT (r - to) In [cjd(t - to)] + 0{£^,e^) , (5.38) 
we are led to the following DRG equation 

^^dfip,t)+e£-VpnF{p)^-2T{t)5f{p,t) , (5.39) 

with r(<) the time dependent relaxation rate, which in the leading logarithmic approximation is given by (for t 3> <o) 

r{t) ^aT [ln{ujDt) + 1] . (5.40) 

The equation for Sf{—p,t) is the same as that above for particles. The DRG equation H5.39|l is the Boltzmann 
equation in the relaxation time approximation. 

// r were a constant the solution to the above equation would be 

1 - e~2rt 

Sfip,t) = -e£-VpnF{p) . (5.41) 

Using the relations given by eqs. ((5.12|) - ((5.13|l with this solution we would then find 

S" (0 ^ -4«7 ^ P V ^ e-' . (5^42) 

This is exactly the result that we would find if we considered the fermion propagator obtained from the Dyson 
resummation with a constant width, as discussed in Sec. IIVBI Thus we see that in the case in which the width is 
constant, the DRG resummation in the relaxation time approximation is equivalent to the Dyson resummation of the 
fermion propagator [23, • The result (|5.42ll leads to a steady-state, drift current at long times given by 



/ T-T^ P £ -VnUFip) = 



J' = -^ I 77;^ f ^ • Vpn^^(p) = ^£\ (5.43) 



and to the electrical conductivity 



a = , (5.44) 

isr ' ^ ^ 

which is in complete agreement with that of the analysis of Sec. lIVBl for constant T. 

However, for T{t) given by eq. H5.40|l the dynamical renormalization group equation (|5.39() leads to a very different 
picture. The solution of eq. H5.39|l with eq. H5.4()|l is given by 

Sf{p,t) = -e£-Vp UFip) e-2"Ttln(^,,t) /"* ^2<,Tt' In(unt') ^^f ^ (5 45) 

Jo 

where we have assumed that the electric field was switched-on at t ~ and that Sf{p,t = 0) = 0. Taking the 
time derivative of this expression, using the identifications eqs. (|5.12|l - (|5.13l) and that the Boltzmann equation for 
antiparticles (with momentum — p) is the same as that for particles (with momentum p), we find to 0{a) 

Q'' it) = 4e2 J ^ f f [1 - 2aTt Hujnt)] = -5,, ^ ^ ^aTt In(^ci)] , (5.46) 
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namely, the solution of the Boltzmann equation in the relaxation time approximation reproduces the leading secular 
term to 0{a) found in the perturbative expansion given by eq. (|4.11|) . 

Thus we reach an important result of this study: the relaxation time approximation in the Boltzmann equation is 
equivalent to only considering self-energy corrections to the internal fermion propagators in the photon polarization. 
This is a reassuring result and establishes a direct link between the perturbative expansion of self-energy corrections 
to the photon polarization and the Boltzmann equation in the relaxation time approximation. 

Of course, these secular terms are an artifact of the approximation, the full solution given by eq. (|5.45|) is bound at 
all times. Its asymptotic long time behavior can be found as follows. First it is convenient to split the time integral in 
eq. as J* dt' = /g dt' + /j*/^^ dt'. In the first integral t' ln(w_Dt') < hence we can replace e2"'r*i"('^«*) - 1, 

leading to an exponentially small contribution at asymptotically long time. The second integral can be computed by 
successive integration by parts which lead to the asymptotic expansion. 



1 



2aTt \n{uJDt) 



(5.47) 



namely the drift current vanishes asymptotically, hence the logarithmic anomalous damping leads to a vanishing 
electrical conductivity. The main result of this analysis is that while a constant damping rate T would lead to a finite 
conductivity given by eq. (|5.44|) , an anomalous logarithmic time dependence of the damping rate leads to a vanishing 
electrical conductivity in the asymptotic, long-time limit. 

Again we highlight that the result given by eq. H5.47|l is a consequence of the relaxation time approximation which 
is equivalent to keeping only self-energy corrections. This result does not apply to the hot QED plasma because 
the vertex corrections introduce the cancellation analyzed in the previous section. The correct expression is given 
below by eqs. H5.59|l and H5.78|) . We derived here eq. H5.47|l to show what would be the solution if the relation time 
approximation were valid, which is not the case here. 



C. Full and linearized Boltzmann equation 

The above discussion on the relaxation time approximation highlights that this approximation is equivalent to 
considering only the self-energy corrections to the photon polarization. This conclusion is based on the equivalence 
of the secular terms between the perturbative result (|4.11|) and the lowest order in a term in the exact solution 
of the Boltzmann equation in the relaxation time approximation (|5.4t)|) . The anomalous logarithms in time are a 
consequence of ultrasoft photon exchange leading to the threshold infrared divergence. 

Just as in the perturbative case we expect that the cancellation of infrared divergences between self-energy and 
vertex corrections will also be manifest in the Boltzmann equation. Thus, we assume for a moment, and will be 
confirmed below, that the coefficient of the cosine term JC{f{., to), /(., tQ):po,p) in eq. (|5.27|) is non-singular and finite 
at po = P and use the result eq. H5.29|l to extract the asymptotic long time limit. In this case the secular term is linear 
in time, and we recognize that the collision term that enters in the Boltzmann equation (|5.33|l is given by 

Cs[f{.,t)J{.,t)]=Ctnfi;t)Ji.,t)] + Cf[fi.,t)Ji.,t)] , (5.48) 

where 

- /■ d^n 

[/(.,i),/(.,i)] = 27re' j ^ J dqo{[K+{p,Ci)*pLiqo,q) + K+ip,q)*pTiqo,q)] Sip - qo - \p + <l\) 

X [[1 - /(P, t)] /(p + q, t) UB (go ) - /(P, i) [1 - /(P + q, t)] [1 + ub (qo)]] 
+ [Kiip,q) *pL{qQ,q) + K^{p,q) V(9o,g)] S{p - qo + \p + q\) 

X [[1 - /(p, t)] [1 - 7(-p - q, t)]nBiqo) ~ /(p, i)7(-P - q, t)[l + ns(go)]] } , (5.49) 

— f d^n 1 1"^°° 

Cf[f{.,t)J{.,t)] = Tre^ / \ _ / dgo{[i^3+(p,q) >+(go, g) + i^3+(P> q) V('?o,'z)] 5{p - qo - \p - c^\) 

J i^'^'j IP qi j-oc 

X [[1 - /(p, t)] npiqa) nsdp - q|) - /(p, t) [1 - npiqo)] [1 + n^dp - q|)]] + 
+ [K^{v.Oi) *P+iqo,q) + K;^{p,q) *p^{qo,q)] S{p + qo + \p - q\) 

X [[1 - /(p,t)] [1 - npiqo)] [1 + ubHp - q|)] - /(p,t) npiqo) ubHp - q|)] } . (5.50) 

The reason that in the soft-fermion contribution there appear the free fermion distribution functions np{qQ) and 
not / is a consequence of the fact that the HTL correction to the vertex does not lead to secular terms to lowest order 
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as discussed above and analyzed in detail in Ref . [TP| . This will become clear below when we identify self-energy and 
vertex contributions in the collision kernel of the Boltzmann equation and we compare the perturbative solution of 
the Boltzmann equation with the perturbative expansion of the function Q{t) studied in the previous sections. 

We notice that the delta functions S{p — go + |p + q|) in eq. (|5.49l) . which multiplies terms that mix particles (/) 
with antiparticles (/), and d{p + qo + \p + q|) in eq. H5.50|l can only be satisfied ii qo ^ T since p ^ T and q ~ eT. 
However, in obtaining the kinetic current we have neglected precisely these type of processes which are responsible 
for very fast oscillations on time scales ~ 0{1/T) [see discussion after eqs. (|5.3ll and (|5.4|) ]. 

To be consistent with the kinetic description, we must consider only processes in which the transfer of energy 
and momentum is go, g ^ T and hence neglect the second terms in eqs. (|5.49|l and (|5.50|l . Therefore, the correct 
Boltzmann equation to this order is given by 

-/(p,t) + e£-Vp/(p,i) - 27re2y ^ J dqol[K+{p,q) *pL{qo,q) + K+{p,q) *pT{qo,q)] 

X [[l-./(p,<)] ./(p + q,<) nBiqo)- f{p,t) [1 - /(p + q, i)] [1 + neiqa)]] 
X 6{p -qo-\p + q|) + I \ I [K^ip, q) *P+{qo, q) + K;^{p, q) *P-{qQ, q)] 

^IP qi 

X [[1 - /(p, t)] npiqo) ubHp- q|) - /(p, t) [1 - npiqo)] [1 + n^dp - q|)]] 

X S{p-qo-\p-q\)} . (5.51) 



The Boltzmann equation for antiparticles (with momentum — p) is exactly the same as eq. (|5.51|l but with the 
replacement f{p,t) — > /(— p,i) on both sides of the equation. 

The linearized approximation of the Boltzmann equation H5.51|l can be obtained by setting 

f{p,t)^np{p)+df{p,t) , J(-p^t)^nFip)+67i-p,t) , (5.52) 

where Sf and Sf are the departure from equilibrium, and keeping only terms linear in Sf and 6f in the collision 
term. Since the collision term vanishes for the equilibrium distribution (np), we find the final form of the linearized 
Boltzmann equation 



d_ 
di 



6fip,t) + ee-Vpnp{p) = -S^e^y ^ dqo!^[K+ip,q) *pL{qo,q) + K+{p,q) *pT{qo,q)] 

X [(5/(p, t) [1 + nsiqa) - npUp + q|)] - Sf{p + q, t) [ns(go) + npip)]] 
X d{p-qo - |p + q|) + —^—[K^ {p,(i) *p+{qo, q) + {p,q) *P-{qo,q)] 

^IP qi 

X Sf{p, t) [1 + nsdp - q|) - npiqo)] S{p - go - |p - q|)l . (5.53) 



The departure from equilibrium for antiparticles Sf{—p,t) satisfies a similar equation but with the replacement 
^f{Pit) ^f{~Pit)- We note that the linearized approximation entails that Sf, 6f cx £, hence the replacement 
eS •Vp/(p, i) ^ e£ ■'Vpnp{p) on the left-hand side of the Boltzmann equation is justified because the neglected 
term is of second order in £. 



D. Connection with perturbation theory: the emergence of secular terms 

The linearization of the collision term in the Boltzmann equation has a simple and natural interpretation in terms of 
Feynman diagrams. We recall that the collision term to lowest order is given in Fig.|3 where the fermion propagators 
are those in the external background field. The linearization corresponds to expanding the propagators to linear order 
in the background. There are two types of insertions: (i) The external background field is inserted in the fermion 
lines outside the self-energy (two diagrams) , this corresponds to the relaxation time approximation which is manifest 
by the first term in eq. H5.53|l [proportional to Sf{p,t)] and is depicted in Fig. [rUT al. (ii) The insertion is in the 
fermion propagator inside the self-energy, which corresponds to the second term [proportional to Sf(jp + q, t)] and is 
depicted in Fig. llOf b^. Obviously this last term corresponds to the vertex correction. In fact these diagrams must 
be compared to those of the self-energy and vertex corrections in the perturbative expansion given in Figs. 13 and 01 
respectively, with the external wiggle line in Fig. IIUI corresponding to the external photon line in Figs. |21 and 21 
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(b) 

FIG. 10: Diagrams contributing to the linearized Boltzmann equation. There are two diagrams of type (a) corresponding to 
insertions to each side of the self-energy. The solid line denotes the free fermion propagator without the background gauge 
field. The wiggly line with a circled cross attached to the end denotes to an insertion of the background gauge field. Otherwise, 
the notation is the same as that in Fig. |5| 

In order to see this relation more clearly, we now show that the perturbative expansion of the solution of the 
linearized Boltzmann equation H5.53|l reproduces the perturbative results found in Sec. IIV Dl above. Let us write the 
solution of the linearized Boltzmann equation in the following form 

Sfip,t) = e <5/(i)(p,i) + Sf^'Hp,t) + ■■■ + e^"-! <5/(")(p,i) + • ■ ■ , (5.54) 

then eq. H5.53|l reduces to 

^5/«(p,t) = -e-pn'p{p) , 

: : (5.55) 
|^J/(")(p,t) - 5C[<5/("-i)(p,t)] , 



where n'p{p) ~ dnp{p) / dp and 5C denotes the collision term of the linearized Boltzmann equation (|5.53|) divided by 
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e^. The set of equations (|5.55() can be solved by iteration starting from the lowest order one. Up to order we obtain 
(5/(1) (p,i) - -t£-pn'p{p) , 

x{£ -p n'p{p) [1 + ni3(go) - nF(|p + q|)] - S • pTq n'p(|p + q|) [nsiqa) + npip)]} 
X S{p-qQ - |p + q|) + — [i^;^(p, q) *p+{qo,q) + {p,q) *p^{qn,q)] 

^IP q| 

X £ • p n'p{p) [1 + nsdp - q|) - npiqo)] S{p - go - |p - q|)| ■ (5.56) 

^From these expressions we can now proceed to compute the current 

•r{t)=4e I ^f5f{p,t)^ (5.57) 

where the factor four accounts for spins, particle and antiparticle. It is more convenient to compute dJ'(t)/dt in order 
to compare to the perturbative expression of the function Q^^{t). Inserting eq. H5.56|l into eq. H5.57|l . we find up to 
order that 

/+00 r 
dqo \ iP' - f+q) [Ktiv. q) *PL{qo,q) + K+{p, q) *Pt('?o, ?)] 
-oo L 

X [nB(go) + »F(|p + q|)]+ I _ . [g3+(p,q) *p+{qo,q) + K^{p,q) *p-{qQ,q)] 

^IP qi 

X [1 + ubHp - q|) - npiqa)] \ S{qo - p • q) . (5.58) 



In the above expression we have combined terms of 0{e'^) that arise from the soft- photon contribution. This is done 
by relabelling qo —qg and using the properties 1 + nB{—qo) — —nB{qo) and pL.Ti—qoi q) = —PL,T{qa, q) in the first 
term of ©(e^), as well as by relabelling p —p — q in the second term of this order. Furthermore, we have also used 
q <^ p ^ T to approximate the argument of the delta functions. 

We can now compare the above perturbative expression for dJ'^{t)/dt to that obtained from t/'^ (t) in Sec. IIVDI 
With the identification given by eq. (|2.28() it is clear that the term of O(e^) in eq. (|5.58|l exactly reproduces the 
HTL result eq. (|3.3(l and that the terms of 0{e*) exactly reproduces the two-loop results for the soft-photon and 
soft-fermion contributions given by eq. (|4.28l) . Thus we see that the perturbative solution of the linearized Boltzmann 
equation reproduces the leading secular terms of the perturbative expansion resulting from the self-energy and vertex 
corrections. The cancellation of the contribution of ultrasoft photon exchange observed at the level of the perturbative 
expansion is a direct manifestation of the transport time scale that enters in the Boltzmann equation. This exercise 
not only shows this important result in detail, but it also highlights the following important identifications: 

(i) The term with Sf{p,t) in the collision term of the linearized Boltzmann equation (|5.53l) is associated with the 
self-energy correction, as manifest in the previous study of the relaxation time approximation. 

(ii) The term with Sf{p + q, t) in the collision term is associated with the vertex correction. The cancellation 
between the self-energy and vertex corrections in q — > limit, is manifest in the linearized Boltzmann equation 
1)5. 53(1 as the cancellation between the terms with df{p, t) and 5f{p + q, t) [upon using the fact that the photon 
spectral functions *PL,Tiqo, q) are odd functions of go]- 

This explicit comparison between the perturbative evaluation of the photon polarization and the perturbative expan- 
sion of the solution of the linearized Boltzmann equation (in real time) provides also a justification for ignoring the 
HTL correction for the vertex [see Fig. 0fb)] at the level of the collision term. This diagram does not contribute to 
the leading logarithm for linear secular terms in the perturbative expansion, and hence will not contribute leading 
logarithms to the collision term of the linearized Boltzmann equation. 

We now argue that the equivalence between the perturbative solution of the linearized Boltzmann equation and 
the perturbative evaluation of the photon polarization in quantum field theory in real time can be generalized to all 
orders in perturbation theory, providing a direct identification of the infinite diagrams in the photon polarization that 
are resummed by the linearized Boltzmann equation. This is achieved by the following observation: 
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FIG. 11: Diagram corresponding to the lowest order perturbative solution 5/'^' of the linearized Boltzmann equation. 

(i) The lowest order solution 5f^^^ features linear secular term and corresponds diagrammatically to the HTL-type 
diagram with hard fermion loop, as depicted in Fig. 1111 

(ii) The second order solution Sf'^^\ which is determined by SC[Sf^^^, features quadratic secular term and cor- 
responds to diagrams obtained from the HTL-type diagram in Fig. 1111 by inserting self-energy and vertex 
corrections, respectively, into the hard fermion lines and the vertex attached to the external wiggly line. 

(iii) As discussed in detail in the previous sections, the secular terms with powers i" arise from singular denominators 
of the form (k • p — cj)^" for n = 1,2... in the limit k; w 0. In turn these denominators arise from the pinch of 
two propagators on opposite side of the loop. Thus the iteration above which features powers t, t^, ... arises from 
the insertion of uncrossed ladder type diagrams with self-energy corrected side rails. This is the manifestation 
in real time of the argument provided in references |12L that the Boltzmann equation to leading order sums 
(uncrossed) ladder type diagrams with self-energy corrected side rails. 

^From this observation and the iterative nature of the perturbative solution [see eq. (|5.55() ]. we conclude by induction 
that the n-th order solution 5/^"' features secular term of the form t" and the corresponding diagrams is obtained from 
those of the (n — l)-th order solution i5/^"~^' by performing the same insertion described above. Since the insertion of 
the vertex correction involves only the vertex attached to the external wiggly line, the only set of diagrams that will 
be generated are uncrossed ladder diagrams with rainbow type self-energy insertion in the side rails. Consequently, 
this analysis reveals that the linearized Boltzmann equation H5.53|l resums the uncrossed ladder plus the insertion of 
rainbow type self-energy diagrams in the photon polarization in accord with the results of reference |l8l Il9l | . 

There is an important corollary of the equivalence between the perturbative solution of the linearized Boltzmann 
equation and the perturbative evaluation of the photon polarization in quantum field theory in real time. Namely, 
by identifying the contributions at different orders to the linear secular terms in the perturbative expansion of the 
function G{t), one can understand precisely what are the diagrams that contribute to the collision term in the linearized 
Boltzmann equation. In particular, this equivalence allows a direct identification of the contribution from self-energy 
and vertex corrections to the polarization to the different terms in the linearized Boltzmann equation. Thus, obtaining 
the linear secular terms in Q(t) yields direct information on the collision term of the linearized Boltzmann equation 
to all orders in perturbation theory. 

Furthermore, this discussion clearly highlights that having proved that the contributions from self-energy and vertex 
insertions that leading to secular terms fulfill the Ward identity in section (|III El) , the subsequent series resummed 
by the dynamical renormalization group or Boltzmann equation is fully gauge invariant. 

Another important issue that we want to highlight at this point is that the HTL-resummed spectral functions *pt,l 
and *p± are obtained from the Dyson resummation of the hard thermal loop result, namely these spectral densities 
are function(als) of the hard fermion distribution functions /(p, t) and /(p, t). Thus the question arises: Why did we 
not account for a change in the fermion distribution functions that enter in the HTL-resummed spectral functions? 
The answer to this important question is the following: consider such a variation in the hard fermion distribution 
functions that enter in *pl,t and the resulting contributions are depicted in Fig. 1121 For *pl,t such contribution 
would correspond to the vertex correction with a HTL-resummed three-photon vertex [see Fig. I12f a)]. However, 
such vertex is forbidden by Furry's theorem, a consequence of CPT invariance in a charge neutral medium (i.e., 
vanishing fermionic chemical potential) in the absence of a magnetic field. This conclusion is obviously supported 
by the perturbative expansion, where the absence of such three-photon vertex is manifest. For *p± such contribution 
would correspond to the HTL correction for the vertex [see Fig. I12f bl] which is subleading at the leading logarithmic 
order as we have argued in Sees. HvT^ and FlV Dl Thus to leading logarithmic order the departure from equilibrium 
of the fermion distribution functions are only associated with the hard fermion propagators that generate the kinetic 




FIG. 12: Possible contributions to the linearized Boltzmann equation resulting from a change in the fermion distribution 
functions that enter the HTL-resummed propagators in (a) the soft-photon contribution and (b) the soft-fermion contribution 
(see Figs. El and Cnj. 



current, and not with the internal fermion propagators that enter in the HTL corrections to the soft photon and 
fermion propagators. 

The Boltzmann equation provides a resumniation of these secular divergences. As argued above, by carrying out 
the perturbative expansion of the Boltzmann equation to higher order and identifying secular terms of the form 
with inverse powers of frequency in the Fourier transform, one can identify what type of the diagrams 

in the polarization and which parts of the diagrams (e.g., leading logarithms, etc.) that are being resummed. The 
resummation of diagrams furnished by the Boltzmann equation can thus be put on a similar footing as the resummation 
of diagrams in the usual renormalization group, where by solving the equation for the running coupling via the 
renormalization group beta function and expanding this solution in a naive perturbative expansion, one can recognize 
what type of diagrams and which part of the diagrams are being resummed. Up to the leading logarithmic order 
considered, the linearized Boltzmann equation resums the leading logarithmic contributions from the uncrossed ladder 
plus the insertion of rainbow type self-energy diagrams in the photon polarization [iM [T9l | . 

The resummation of the linear secular term in real time is also akin to computing the decay rate in time dependent 
perturbation theory. By using Fermi's golden rule to extracts the leading linear secular term at long times and expo- 
nentiating the secular term, one obtains the exponential decay of probabilities. This exponentiation is a resummation 
of the perturbative series. 



E. Fokker-Planck equation 

For soft momentum transfer q <C p T we can expand the integrand in the collision term of the linearized 
Boltzmann equation (|5.53l) in powers of q/p, as we did previously in Sec. lIVDl to obtain the leading logarithms in 
the perturbative expansion. This expansion, valid for small momentum transfer is akin to that used in astrophysics 
to obtain the Kompaneets equation from the kinetic equation for the photon distribution function, which is widely 
used to study inverse Compton scattering by hot electrons |3l|. The result is a partial differential equation for the 
distribution function, a Fokker-Planck equation in momentum space. 

The first step in the derivation consists in writing the departure from equilibrium as in eq. H5.12(l in the following 
form 

dfip,t)^n'pip)e-pA{p,t). (5.59) 
The main reason for this definition is that the inhomogeneous term in the Boltzmann equation is proportional to 
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FIG. 13: Diagrams that yield the leading logarithmic contribution to the Boltzmann equation. 

n'p{p) £ ■ p. Our analysis of the leading logarithmic contributions in Sec. lIVDl clearly showed that in order to extract 
the leading logarithms, the transverse and longitudinal photon spectral functions can be replaced by eq. (|4.3U|I and 
the momentum integrals must be carried out with an upper cutoff of order T and a lower cutoff of order eT. These 
spectral functions are the perturbative ones, corresponding to just one hard fermion loop in the soft photon and fermion 
propagators, not the full Dyson series of bubbles. The diagrams that yield the leading logarithmic approximation are 
depicted in Fig. 1131 These diagrams are identified with those of Fig.|Slin the perturbative expansion with the external 
wiggly line in Fig. 1131 corresponding to the external photon line in Fig. 

For g <C p the argument of the delta functions in the collision term in eq. (|5.53|l are approximated by 



|p±q| -P = < ±9P • q 



(p-q)' 



(5.60) 



which simplifies the integral over go in the collision term in eq. H5.53|l . The resulting integrand in the collision term 
can be expanded in powers of q/p. For the soft-photon contribution the leading logarithm arises from terms linear in 

q/p, which leads to a momentum integral of the form jjj, dq/q = ln(l/e). Terms of 0(1) in q/p are odd in p • q thus 
vanish after angular integration, while terms of higher order in q/p lead to higher powers of a but not logarithms. 
For the soft- fermion contribution the leading logarithm arises from terms of 0(1) in q/p, which leads to the same 
momentum integral as for the soft-photon case. Likewise, terms of higher order in q/p lead to higher powers of a but 
not logarithms. 

After some lengthy but straightforward algebra we find the linearized collision term in eq. (|5.53|l to be given by 



C[A] 



irp3 



247r 



In - 5 -pn'pip) 



^"{P,t) 



tanh^j A'(p,<) 



— ^-TT coth — 
4pr2 2T 



A(p,t) 



(5.61) 



where the primes stand for derivatives with respect to p. The Boltzmann equation in this approximation becomes the 
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Fokker-Planck equation 



3y3 



■In - 

247r V e 



It proves convenient to introduce the transport time scale 



A"(p, t) + f - - tanh ^) A'ip, t) ( coth ^ + ^] A(p, t) 



V4pT2 



2r p2 



= -1 



247r 



e4rin(l/e) 



and the dimensionless variables. 



in terms of which the Fokker-Planck equation becomes 

^{x,t) - $"(.T,r) - [ - - tanh I j $'(a;,r) 



t 



P 

X = — 

T 



i^C0th- + -)$(.T,T)=-l 



(5.62) 



(5.63) 



(5.64) 



(5.65) 



where the dot and prime stand for derivatives with respect to t and x, respectively. This equation clearly reveals that 
the time scale for transport is determined by itr defined in eq. H5.63|l . The real-time approach leads directly to this 
Fokker-Planck equation in the leading logarithmic approximation, which reveals at once the transport time scale and 
includes the time dependence of the departure from equilibrium. Thus this approach has a distinct advantage over 
the ladder resummation program advocated in Refs. [m . IT^ [l9j which only captures the asymptotic, steady-state 
solution. 



F. Solution of the Fokker-Planck equation: asymptotics 

The solution of the Fokker-Planck equation H5.65|l can be written as 

$(a;,r) = (I){x,t) +<^oo{x) , 
where <f>{x,T) is a solution of the homogeneous time-dependent Fokker-Planck equation 

2 



{x,t) - (j) {x,t) 



2 x\ / 3 X 

_ _ tanh - j 0'(a:, r) + ( — coth - + -2 



x,t) = , 



^x 2 J \Ax 

while, the steady-state solution $oo(a;) obeys the time-independent inhomogeneous equation, 



*'c^(2;)+ 0-tanh0 $;^(a;)- (^Acoth| + A^$^(a;) = i 



(5.66) 



(5.67) 



(5.68) 



Eq. (|5.67|l is a parabolic equation that can be solved by expanding 0(a;, t) in terms of the eigenfunctions ^„ (x) satisfying 















— tanh — 




dx'^ 




2y 


dx 



in the interval < x < c» with ?/'n(0) — ipnioo) = 0. Setting 



i^n{x) = COsh^ , 

X 2 



we can rewrite eq. (|5.69|) in a Schrodinger-like form 

d^ 



dx'^ 



+ V{x) 



£,n{x) — \n£,n{x) 



where 



T^/ N 1 1 1 , X 3 , I 2 

V(x) = 5 tanh - + — coth - + ^ 

4 2 cosh^ fx 2 4x 2 x2 



(5.69) 



(5.70) 



(5.71) 



(5.72) 
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This potential has a repulsive singularity for x — > 0, 



V{x) 



2x2 



and for x 3> 1 approaches its asymptotic value with an attractive Coulomb tail 



V{x) has its minimum at x = Xyam — 15.997 . . . with 

V"(x,ni„) = 0.2421873 . 



y"(xmin) = 0.0000609. 



(5.73) 



(5.74) 



(5.75) 



In addition, ^(x) > ^(xmin) and V"{x) > for < x < oo. 

The Coulomb-like tail of the potential suggests that the spectrum of eq. H5.69|l is formed by an infinite number of 
discrete eigenvalues with an accumulation point at A — 1/4 — and a continuum that starts at A = 1/4. Therefore 
the discrete eigenvalues Afe (with A; = l,2,3,...) are in the interval V^(xinin) < < \- Thus the solution of eq. (|5.()7|) 
can be written as, 



r) = ^ Cfe e-^" ^ Mx) + dX c(A) e-^ ^ Va(^) 



(5.76) 



where the Fourier coefficients Ck and c(A) are obtained from the initial data (/)(x, 0) by scalar product with the 
eigenfunctions. The late times asymptotics is dominated by the ground state. 



(x,t) 



ipiix) 



(5.77) 



A simple estimate from eq. (15 .751) . namely, Ai « V"(xniin) + \/V" {xmin)/'^ suggests that Ai « 0.247664. Therefore, 
A{p,t) approaches the steady-state solution as 



Aip,t) A(p,cx3) 



24 7rci 



-t/(4.038... ttr) 



e3rin(l/e) 

The steady-state drift current at asymptotically large time t 3> itr is given by 



ipiip) ■ 



n'p{p) A{p) 



(5.78) 



(5.79) 



where the factor 4 accounts for spin and the fact that the departure from equilibrium for the particle and antiparticle 
distribution functions are the same. In terms of the steady-state solution ^oo{x) and the dimensionless quantities 
introduced above in eq. H5.64|l . the conductivity is finally given by 



16T 



Tre^ ln(l/e) Jq 



dx x^ n'p{x) ^ooix) ■ 



^From the differential equation H5.68|l we can write the following identity 



*'4(^)+*'oo(^) 



2 , X 

tanh — 

X 2 



3 , X 
— coth — 
4x 2 



(5.80) 



(5.81) 



which is now inserted in the expression for the conductivity eq. H5.80|l . After integration by parts we finally find that 



32 T 



ne^ ln(l/e) Jq 



dx n'p{x) 



l + ^coth| ) $^(x)-f x2 ^^{x) 



(5.82) 



It is straightforward to confirm that this expression is exactly 2(5[x]/3 given by (4.3) in Ref. [3 after the rescaling of 
the function x{p) iii that reference to our dimensionless function <^ryo{x). Considering the conductivity as a functional 
of $00(2;) given by eq. (|5.82l) . we find that the variational condition 



(5.83) 
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leads to the steady-state Fokker-Planck equation (|5.68() . 

Clearly an analytic solution of the steady-state Fokker-Planck equation H5.ti8|l is not available. However, before we 
embark on a numerical study it is important to understand the small and large x behavior of the solution. For this 
purpose we use the small and large x expansion of the hyperbolic functions. After some algebra we find that the 
solution regular at the origin has the power expansion 



H<{x) 
U{x) 



a iJ< (x) + /< (x) , 

x" [1 + ai x^ +a2 x'^ + 0{x^)] 



110 



157 
254100 



x^ + Oix'^) 



(5.84) 



ai^ ^':i^ = 0.0717902... 



where 

For a; ^ 1 the solution that does not grow exponentially is given by 

$>(x) =&ff>(a;)+/>(a;) 



6099- 1526 V15 
= = 0.00223517. 



(5.85) 



(5.86) 



where Hy(x) has an asymptotic expansion in positive powers of 1/x and e ^. Neglecting the exponentially small 
corrections, we find 



H^{x) x^[\A:X^ [l + 0(e-")] 

i-fV^^r:^ 

k>l \ / k \ / k 



35 385 
^+16^+512^ + ^(^ ) 

Here, 4'(a, 6; z) stands for a confluent hypergeometric function 33] and (/i)fc 
I>{x) has the asymptotic behavior 



(5.87) 



/>(.t) = -jX 



+ 1) ■ ■ ■ {fi + k — 1) . The function 



(5.88) 



The functions H^{x) are the corresponding solutions of the homogeneous equation, while I^{x) are the particular 
solutions of the inhomogeneous equation that are finite at the origin and do not grow exponentially at infinity. The 
asymptotic linear growth with x of J> (x) is a direct consequence of the inhomogeneity. In addition, the homogeneous 

equation has a linearly independent (irregular) solution that grows as x ^ for x 0. 

The coefficients a and b are found by integrating both the homogeneous and inhomogeneous equation with the 
initial conditions gleaned from eq. (|5.84|) on the function and its derivative at the origin for i?< (x) and /< (x) and 
matching the function and derivative to the asymptotic form eqs. (|5.87ll - (|5.88ll at a value of cc ^ 1. Such procedure 
is numerically straightforward, we find that the asymptotic form is achieved for relatively small values of a; ~ 3 and 
the coefficients are found to be given by a = —0.896 . . ., b ~ —0.01 . . .. The numerical solution is shown in Fig. ITTI 
With this solution we finally find the following value of the conductivity 



15.698 T 

e2 ln(l/e) ' 



(5.89) 



which is consistent with the numerical value quoted in Ref. 9], a — l^fn^i}!) ■ However we point out that the 
asymptotic forms for the solution for small or large x are very different from the variational basis chosen in Ref. , 
in particular near the origin the solution is not an analytic function of x. 



G. Beyond leading logarithms and the Boltzmann equation 



In Sec. IIVDI we have discussed the diagrams that will contribute to next to leading order, which are necessary to 
fix the constant inside the logarithm for the transport relaxation time and the conductivity. These diagrams will yield 
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FIG. 14: Numerical solution of the steady-state Fokker-PIanck equation 15.681 . 



new contributions in the collision term. To next to leading order the nonequilibrium fermion distribution function / 
will replace the equilibrium one np in the three-loop diagram for the polarization in the soft-fermion contribution ^ 
[see eq. (|5.5U|l ]. the collision kernel must also include the first diagram in Fig.|Sl the second diagram will arise from 
linearizing in the departure from equilibrium as discussed above. 

The discussion after eq. 1)5.4(1 indicates that the validity of the Boltzmann equation for the evolution of the particle 
and antiparticle distribution functions is reliable provided processes in which particles and antiparticles mix can be 
neglected, namely, processes involve momentum exchange q <^ T. However, if there are processes that mix particles 
and antiparticles that lead to secular terms in the perturbative expansion, then the contribution from the spin current 
in eq. ((5.4(1 must be accounted for in the evolution equation. In this case one would have to find the equation of 
motion for the bilinears these equations will likely involve the correlators (6^5) leading to a coupled system 

of equations. The new terms are associated with spin precession and in general must be included for a consistent 
nonequilibrium description '32*1. The perturbative analysis and the direct correspondence of the secular terms from 
the perturbative solution of the Boltzmann equation indicate that such new contributions do not appear at leading 
logarithmic order. Thus the simpler Boltzmann equation for the particle and antiparticle distribution functions in 
this leading logarithmic approximation can be interpreted as a low-energy limit of the full set of evolution equations, 
in the sense that processes that mix particles and antiparticles are neglected. 

VI. CONCLUSIONS 

In this article we began a program to study transport phenomena in ultrarelativistic plasmas via the dynamical 
renormalization group resummation. We focused on the DC electrical conductivity in a ultrarelativistic QED plasma 
and extracted the conductivity from the long time behavior or a kernel directly related to the retarded photon 
polarization. 

We began by studying the long-time behavior of this kernel in perturbation theory in real time. Pinch singularities 
in the limit of long wavelength and low frequency that are ubiquitous in the usual perturbative expansion of the 
imaginary part in the polarization in frequency-momentum space, are manifest as secular terms in real time, namely 
terms that diverge in the long-time limit. At any finite time the perturbative contributions to the kernel are finite but 
the long time behavior of the perturbative expansion is divergent reflecting these pinch singularities. We first extracted 
the leading secular terms at lowest order in the perturbative expansion which includes both self-energy and vertex 
corrections for the fulfillment of the Ward identity. While the self-energy correction features a logarithmic singularity 
on the fermion mass shell arising from the exchange of ultrasoft photon which leads to the anomalous damping rate 
of the hard fcrmions in the loop, the addition of the vertex correction, as required to fulfill the Ward identity cancels 
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the contributions from the anomalous fermion damping. This cancellation is at the heart of the difference between 
the quasiparticle and transport relaxation time. We then extracted the leading logarithmic contributions to the 
leading secular terms and identified the diagrams that yield the leading logarithms and those that yield subleading 
contributions without logarithms of the coupling. 

After the study of the secular divergences in the pcrturbative expansion, we introduce a resummation program 
via the dynamical renormalization group applied to the equations of motion of the gauge invariant single-particle 
distribution functions of hard fermions. The resulting dynamical renormalization group equation in real time is 
the Boltzmann equation. By expanding the solution of the linearized Boltzmann equation in perturbation theory, we 
established a direct link between the perturbative expansion of the polarization and the linearized Boltzmann equation 
by explicitly showing how the Boltzmann equation reproduces the secular terms in a perturbative expansion. This 
comparison allowed to identify the different terms in the linearized Boltzmann equation with the self-energy and 
vertex corrections to the polarization. This detailed study leads to a deeper understanding of the resummation 
of the perturbative series in quantum field theory via the Boltzmann equation. Furthermore, it allows to identify 
unambiguously the diagrams that contribute to the collision term of the linearized Boltzmann equation. Just as the 
renormalization group equation in deep inelastic scattering or critical phenomena provides a resummation of parts 
(such as leading logarithms) of select type of diagrams, the Boltzmann equation also resums parts (leading logarithms) 
of select diagrams (uncrossed ladder plus rainbow type self-energy corrections). 

Thus the dynamical renormalization group provides directly in real time a link between the quantum field theoretical 
and the Boltzmann kinetic approaches to transport phenomena. Furthermore, we have explicitly shown how the 
relaxation time approximation of the Boltzmann equation is equivalent to neglecting the vertex corrections in the 
linearized approximation and features the anomalous logarithmic time dependence manifest in the hard fermion 
damping rate. The conductivity vanishes in the relaxation time approximation as a consequence of the anomalous 
logarithms in time. As discussed in Sees . IIV Bl and IVBl such an approximation is not valid for the hot QED plasma. 

Recognizing that the leading logarithmic contribution to the collision term arises from the region of momentum 
exchange eT < q < p T, we expanded the collision kernel of the linearized Boltzmann equation in powers of q/p to 
extract the leading logarithmic contribution. We then obtained a Fokker-Planck equation in momentum space for the 
small departure from equilibrium. This equation clearly reveals that the transport time scale is itr = e*T^\n{i/e) ■ 
time dependent Fokker-Planck equation is solved by expanding in the eigenfuntions of a positive definite Hamiltonian. 
For late times, its solution approaches the steady-state solution as ~ g-t/(4.038...ttr)^ solved analytically the 
steady-state Fokker-Planck equation for small and large momenta. We matched numerically the solutions and found 
that the DC conductivity is given by cr = to leading logarithmic order. This result is within less than 0.1 % 

that quoted in Ref. 9]. We established contact with the variational formulation of Refs. [ol. ITsL ITol by showing that 
the steady-state Fokker-Planck equation can be obtained from a variational functional which is the DC conductivity. 
We have identified the diagrams that contribute beyond the leading logarithmic order and provided a discussion of 
the corrections that are necessary to be included in order to go beyond the Boltzmann equation. 

The results of this study highlight that the dynamical renormalization group is the concept that provides a direct 
bridge between quantum field theory and the kinetic description in real time and puts transport phenomena on a 
similar footing as critical phenomena, namely, the Boltzmann equation is the renormalization group equation that 
provides a nonperturbative resummation of diagrams. It establishes this formulation as an alternative to studying 
transport phenomena in relativistic plasmas firmly based on the microscopic quantum field theory. In the next step of 
our program we will apply these ideas to study transport phenomena in QCD relevant to the physics of the quark-gluon 
plasma. 
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APPENDIX A: IMAGINARY-TIME PROPAGATORS 



In this appendix we review the imaginary-time propagators in their spectral representations. The fermion propagator 
is given by [g 



S{iuj„,p) = J 



, Pf(po,p) 
dpQ 

OO 



Po - 

where ujn = (2rt + 1)ttT. The spectral function pf{po, p) for hard (or free) fermions is given by 

Pf(po,p) = 7+(p) P+{Po,p) +7-(p) P-{Po,p) , 
7±(P) = (7° T7 ■ P)/2 , P±iPo,p) = S{Po Tp) , 



(Al) 



(A2) 



where p = |p|, p = p/p and p±{po,p) are the free particle (+) and antiparticle (— ) spectral functions, respectively. 
The soft (or HTL-resummed) fermion propagator *S{iuj,p) has the same form as that given in eq. ljAl|l but in terms 
of the HTL-resummed fermion spectral function, 



*Pf(po,p) 7+(p) *P+{Po,p) +7-(p) *P-iPo,p) , 

*P±{po,p) = Z±{p) 6{po - w±(p)) + Z^{p) 5{po + uj^{p)) + (3±{pa,p) Q{p^ ~ pi) , 
where Z±{p) are the residues of the timclikc fcrmionic quasiparticlc poles at ui±{p), respectively, and 



/3±(po,p) 



(A3) 



(A4) 



with = e^T^/S the thermal for mion mass. 

The longitudinal and transverse photon propagators are given by 



DT{iVn,q) = J 



1 



dqo 



PL{qo,q) 
dqo : 

-00 90 - ll^n 

PT{qo,q) 
qa - Wn 



(A5) 



respectively, where Vn = 2'nnT. The frequency independent term is Dl arises from the instantaneous Coulomb 
interaction and is neglected in our calculation because it does not contribute to the imaginary part. For hard photons 
we have = and 



PT{qo,q) = -^l^ilo - 9) - S{qo + q)] 



(A6) 



The soft photon propagators *DL{iVn,q) and *DT[ii'n, q) have the same form as those given by eq. (IA5|I but in terms 
of the HTL-resummed spectral functions 



PL(T){qo,q) = sgn((7o) ZL(T){q) S[qo - uj^q)] +PL{T){qo,q) 0(g^ - ql) 

where 2'i(7')(g) are the residues of the timelike quasiparticlc poles at u>L(T){q) and 



(A7) 



l3L{qo,q) 
f3T{qo,q) 



Imng^^(go,9) 



TT [q^ + RenHTL(<Zo, q)? + [lmUr^{qo, q)]^ 
1 lmW^^^{qo,q) 



IT [ql - g2 _ RenHTL(g„^ ^)]2 + [ininHTL(<zo, g)]^ 



(A8) 



In the above expressions, H^™ (go, q) is the longitudinal (transverse) photon self-energy in the HTL approximation | 



nr^(go,g) 
nr^(go,g) 



6 
12 



2-^ln 



290 , 90 



qo + q 



qo-q 



in 



e'T^qo 
6 q 



qo 






qo 


+ q 




— iiT 


q 






qa 


- q 







0(9' - qo) , 

12 q 



q'o 



1-^ eiq'-ql), 



(A9) 
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The frequency independent part in Di(w„,g) and *DL{ii'n,q) arises from the instantaneous Coulomb interaction. 
Furthermore, we wiU adopt the convention that bosonic Matsubara frequencies are denoted by = 2m7rT and the 
fermionic ones by Um — (2m + 1)ttT. 

The spectral representation of the fermion and photon propagators facilitate the sum over internal Matsubara 
frequencies. The convenient formulas to carry out these sums systematically are given by jSjJ 

T^F(iz/„) = -Y,^esF{z,)nB{z^) , 

TY,P{i^n) = ^Resi^(zO "f(2») , (AlO) 

where Res refers to the residues at the (simple) poles of the function F(z), is the Bose distribution and np is 
the Fermi distribution (for vanishing chemical potential) 

^^M^ eWT_l ' "J^H^ e^/T + i ■ (All) 



APPENDIX B: REAL-TIME PROPAGATORS 



In this appendix we summarize the real-time propagators which are used in the main text to derive the quantum 
Boltzmann equation. The free fermion propagators (with zero chemical potential) are defined by 

(*-(x,t)*^(x',t')) / ^S;''{t,t')e'P<---"> , 

s++it, t') = s>{t, t') 9{t - 1') + s<{t, t') e{t' - 1), 
Sp-it, t') = s>{t, t') e{t' -t) + s<{t, t') 9{t - 1') , 

S^^{tX) = S^{t,t') , (Bl) 

where (■ ■ •) denotes expectation value with respect to the initial density matrix and a, 6 = ± refer to fields in the 
forward (+) and backward (— ) time branches. For hard fermions the Wightman functions read 

5>(t,0 = -*{7+(p)[l-/(p,M]e-*''^*"*n 7-(p)7(-p,to)e*''(*-*')}, 

S<{t,t') - z{7+(p) /(p,io) e-'P(*-*') + 7-(p) [l-7(-p,io)] e^^^'-''^] , (B2) 

where /(p,to) and /(p,to) are the nonequilibrium distribution functions for the hard fermions and antifermions at 

the initial time to, respectively. For soft fermions the HTL-resummed Wightman functions *Sp-{t, t') are expressed in 
terms of the spectral functions as 

-t-CXD 



S>{t,t') = \ dgoXPo,p)[l-";^(po)]e-^^'° 



+ 00 



*S<{t,t') ^ I dqo*pF{po,p)npipo)e''P°^'-''\ (B3) 



where *pf{po,p) is given by eq. I|A3I) . 

The longitudinal photon propagators are given by 



(Ag(x,t)A[;(x',t')) / j;J^Dt^^{t,t') e^'i-(--') , 



(27r)3 

D++ (t, t') = ^Sit- t') + Dl^{t, t') Bit t') + D< ^{t, t') 9{t' - t) , 
Dl~ it, t') = -^ 5{t - t') + Dl^{t, t') 9{t' ^t) + D<^it, t') 0{t - t') , 

Dtl{t,t') = Dl^{t,t'). (B4) 
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The Wightman functions are expressed in terms of the spectral functions as 

/+00 
dqoPL{qo,q) [l + "B(go)] e-^9o(t-0 ^ 
-OO 

/+00 
d9oPL(go,9)ns(go)e-'*(*-*') . (B5) 
-OO 

where pL = for the hard photon propagator and pL = *Pl for the soft photon propagator. 
The transverse photon propagators are given by 

D++{t,t') = D>^{t,t') 0{t - t') + D<^{t,t') e{t' - t) , 
D^-{t, t') = D> Jt, t') eit' -t)+ D< Jt, t') 9{t - t') , 

D$^^{t,t')=Dfjt,t'), (B6) 

where (q) — 5'-^ — q^q^ is the transverse projector. The Wightman functions are expressed in terms of the spectral 
functions as 

/ + 00 
dqoPTiqo,q) [I + riBiqo)] e'^'^^*-*'), 
-OO 

/ + OC 
dqo priqo, q) ns(go) e-**(*-*') , (B7) 
-OO 

where pr = PB for the hard photon propagator and pr = *PT for the soft photon propagator. 
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